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. After the initial discovery of the so-cahed "spin crisis in the parton model" in the 1980's, a large 

5^ \ set of polarization data in deep inelastic lepton-nucleon scattering was collected at labs like SLAC, 

DESY and CERN. More recently, new high precision data at large x and in the resonance region 
have come from experiments at Jefferson Lab. These data, in combination with the earlier ones, 
allow us to study in detail the polarized parton densities, the dependence of various moments 
of spin structure functions, the duality between deep inelastic and resonance data, and the nucleon 
structure in the valence quark region. Together with complementary data from HERMES, RHIC 
and COMPASS, we can put new limits on the flavor decomposition and the gluon contribution to 
the nucleon spin. In this report, we provide an overview of our present knowledge of the nucleon 
spin structure and give an outlook on future experiments. We focus in particular on the spin 
structure functions gi and 52 of the nucleon and their moments. 
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1 Introduction 

The measurement of a spin dependent observable is generally a daunting task, but with rich rewards, 
because spin seems to have a scalpel-like ability to expose weaknesses and failures of theories. Witness 
the switch from S^T io V — A, which led to the Weinberg-Salam electroweak theory, or the demise of 
Regge poles, which had successfully described hadronic cross-sections and shrinking diffraction peaks, 
or the spin crisis in the parton model of deep inelastic scattering, supposedly resolved by an anomalous 
gluon effect, now shown to be untenable, and, most recently, the realization that some of the information, 
gathered for 40 years on the fundamental electromagnetic form factors of the nucleon, is unreliable. 

In this review we attempt to survey the tremendous experimental and theoretical effort, mainly 
involving studies of polarized lepton-hadron scattering and polarized proton-proton reactions, that has 
led to our present knowledge of the internal spin structure of the nucleon. 

Deep inelastic lepton-hadron scattering (DIS) has played a seminal role in the development of our 
present understanding of the sub-structure of elementary particles. The discovery of Bjorken scaling in 
the late nineteen-sixties [1] provided the critical impetus for the idea that elementary particles contain 
point-like constituents and for the subsequent invention of the Parton Model. DIS continued to play 
an essential role in the long period of consolidation that followed, in the gradual linking of partons 
and quarks, in the discovery of the existence of missing constituents, later identified as gluons, and 
in the wonderful confluence of all the different parts of the picture into a coherent dynamical theory 
of quarks and gluons - Quantum Chromo dynamics (QCD). Polarized DIS, involving the collision of 
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a longitudinally polarized lepton beam with a longitudinally or transversely polarized target, provides 
complementary and equally important insight into the structure of the nucleon. 

At first sight the theoretical treatment of the polarized case seems to mimic the unpolarized case, 
with the structure functions Fi^2{x) replaced by the polarized structure functions gi^2{x), and with 
parton densities q{x) replaced by polarized densities Aq{x). But it turns out that the polarized case is 
much more subtle: there is an anomalous gluon contribution to gi{x), and g2{x) has no interpretation 
at all in purely partonic language. 

The latter insights were mainly inspired by the unexpected results of the European Muon Collab- 
oration (EMC) measurement of gi{x) in 1988 [2]. The first excitement caused by this experiment was 
its indication of the failure of a sum rule due to Ellis and Jaffe based on the assumption that the 
contribution from strange quarks to gi is negligible [5]. It was soon realized, however, that there were 
more profound ramifications, which led to an intense scrutiny of the theory, since they implied a "spin 
crisis in the parton model" [1] — the spins of the quarks seemed to provide only a tiny fraction of the 
spin of the nucleon, in contrast to the situation in simple-minded constituent quark models of hadrons, 
where the quark spins account for a very large fraction of the proton spin. The crisis was believed to 
be resolved via a large polarized gluon contribution. 

Of course the parton model predates QCD. In the more general field-theoretic framework we know 
that Bjorken Scaling |5], i.e. the fact that structure functions and parton densities depend only on x, 
cannot hold exactly, and these functions have a Q^-dependence which can be calculated perturbatively 
in QCD, resulting in some of the most stringent tests of the validity of the theory. Moreover, because of 
the unfortunate need to renormalize the theory, the parton densities lose their simple physical meaning, 
and their actual functional form depends upon the renormalization scheme employed. 

At present the situation, as will be discussed, is full of interest. 

• Measurements of the polarized gluon density suggest that it is much too small to resolve the spin 
crisis [HI E]. This almost certainly implies that the partons possess orbital angular momentum, 
and it appears possible to estimate this, at least for the quarks, via a study of deeply virtual 
Compton scattering on protons [8]. 

• More precise data expected from the COMPASS experiment at CERN and the pp program at 
RHIC will allow further scrutiny of the validity of the above conclusion. 

• There are now significant measurements of g2{x) which can be used to test the Wandzura-Wilczek 
approximation [5], the Burkhardt-Cottingham sum rule [TU], and the Efremov-Leader-Teryaev 
sum rule [TT] . 

• Measurements at Jefferson Laboratory are probing the hitherto inaccessible region of large x, 
where there are intriguing predictions about the behavior of the ratios Aq{x)/q{x), as well as 
the region of low Q^, where higher twist effects are important and where the issue of "duality" 
between the resonance and deep inelastic regions can be studied. 

• By combining DIS data with the growing reservoir of data on semi-inclusive DIS (SIDIS) it should 
become possible to learn about the polarized sea densities Am and Ad and to resolve the present 
disagreement between DIS and SIDIS about the sign of the strange quark density As{x) +As{x). 

The aim of this review is phenomenological, i.e., it tries to strike a reasonable balance between theory 
and experiment. The theoretical treatment is thus conventional QCD and does not address interesting, 
and sometimes profound, matters like the Chern-Simons current and non-perturbative effects such 
as instantons, axial ghosts, and the U(l) problem. These are discussed in detail in the the review of 
Anselmino, Efremov and Leader [12], and in the review of Bass [13], who also examines the consequences 
of a fixed pole in the virtual-photon Compton amplitude. 
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1.1 Definitions and Formalism 



Consider the inelastic inclusive scattering of polarized leptons on polarized nucleons. We denote by 
m the lepton mass, k {k') the initial (final) lepton four-momentum and s [s') its covariant spin four- 
vector, such that s ■ = (s' • fc' = 0) and s ■ s = —1 [s' ■ s' = —1); the nucleon mass is M and 
the nucleon four-momentum and spin four-vector are, respectively, P and S. Assuming, as is usually 
done, one photon exchange (see Fig. [1]), the differential cross-section for detecting the final polarized 
lepton in the solid angle dVL and in the final energy range (£", E' + dE') in the laboratory frame, 
P = (M, 0), k = {E, fc), k' = {E', k'), can be written as 



cr 



where q 



dQdE' 

k — k' and a is the fine structure constant 



E' 



lepton 



nucleon 




Figure 1: Feynman diagram for deep inelastic lepton-hadron scattering. 



The leptonic tensor L^y (summed over the unobserved final lepton spin) is given by 

L^u{k, s] k) = 

J2 Hk'^ s') 7^ uik, s)Y [uik\ s') 7. uik, s)] (2) 

s' 

and can be split into symmetric (S) and antisymmetric (A) parts under /i, u interchange: 

V(fc, s; k') = 2{Ll'J {k; k') + iL^^J {k, s; k')} (3) 

where 

LjfJ {k; k') = kX + - 9,u ik-k'- m') 
L^^J {k,s;k') = me^y^fs q". (4) 

The unknown hadronic tensor W^i, describes the interaction between the virtual photon and the nucleon 
and depends upon four scalar structure functions, the unpolarized functions Fi 2 and the spin-dependent 
functions gi^2 (ignoring parity- violating interactions). These must be measured and can then be studied 
in theoretical models, in our case in the QCD-modified parton model. They can only be functions of 
the scalars and q ■ P. Usually people work with 

g2 = _g2 and XBj = Q^/2q ■ P = Q'^/2Mu (5) 

where 1/ = E — E' is the energy of the virtual photon in the Lab frame. Xbj is known as "x-Bjorken", 
and we shall simply write it as x. We also refer to the invariant mass of the (unobserved) final state, 
W = ^{P + qy = ^M^ + 2Miy - Q^. 



4 



Analogous to Eq. Q one has 

W,,iq; P, S) = Wl,^\q; P) + i W^^^Kq-. P, S). 
The symmetric part, relevant to unpolarized DIS, is given by 



g2 St''' 



Mv 



P., 



P.q 



J^u n qp 
g2 



The antisymmetric part, relevant for polarized DIS, can be written as 



(6) 



(7) 



Wl,^\q- P,s) = 2 e^,^p q'^l M^S^diu, Q^) + [MuSf - {S ■ q) ^2(1^, 



(8) 



or, in terms of the scaling functions gi 2 



Wl,^\q;P,s) 



2M 
P-q 



{S ■ q) P^ 

{P-q) 



92ix,Q^) 



Note that these expressions are electromagnetic gauge- invariant: 



q^W. 



0. 



In the Bjorken limit, or Deep Inelastic regime. 



q' = Q' 



00 



u = E-E' 



00 



X fixed 



(9) 
(10) 

(11) 
(12) 



the structure functions Fi^2 and gi^2 are known to approximately scale, i.e., vary very slowly with at 
fixed a; - in the simple parton model they scale exactly; in QCD their evolution can be calculated 
perturbatively. 

Differences of cross-sections with opposite target spins are given by 



7(^' ^' P^ k') - (k, s, P, S; k') 



dn dE 



2Mg4 E 



(13) 



After some algebra (for a detailed derivation see [12]), one obtains from Eqs. ( fTSlfTOi) expressions for 
the following polarized cross-section differences (note that here, and in the following, we have included 
the lepton mass terms that are usually ignored) : 

• For the lepton and target nucleon polarized longitudinally, i.e. along or opposite to the direction 
of the lepton beam, the cross-section difference under reversal of the nucleon's spin direction 
(indicated by the double arrow) is given by 



d^a 



d^a 



dx dy dx dy 



y y'^{M'^x'^ + m'^)\ 2MV?/ 

^"2 ) ^'-—Q^^' 



(14) 
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For nucleons polarized transversely in the scattering plane, one finds 



dx dy dx dy 



16a2 /2Mx\ 



y 



M^x^y^ 



1'^ 



2m?'y 



91 + 92 



Here we have used 



V P ■ q 



(15) 



(16) 



E P-k 

In principle, these two independent observables allow measurement of both gi and g2 (as has been 
done at SLAG and in Jefferson Lab's Halls A and C), but the transverse cross section difference is 
generally smaller because of kinematic factors and therefore more difficult to measure. Only in the past 
few years has it been possible to gather precise information on g2, which turns out to be usually smaller 
than gi in DIS. 

Since experimental results are often presented in the form of asymmetries, which are the ratios of 
cross-section differences to the unpolarized cross- sect ion, we comment here briefly on the latter. The 
unpolarized cross-section is given by 



dx dy xy Q"^ 



xy 



1 - 



2m' 



M^x'^y'^ 



(17) 



It is a non-trivial task to obtain separate information on Fi and F2 and usually data are presented 
for F2 and R where 

' F2 



R=[l+Y 



and 



7 



2xFi 
AM^x'^ 



1 



(19) 



For a longitudinally polarized target the measured asymmetry is 



All 



da — da 



2 da unpaid 

and for a transversely polarized target the measured asymmetry is 

da^^ - da^^ 

A\ = 



2da unpaid 



It is customary to introduce the (virtual) photon-nucleon asymmetries ^1,2, 



(20) 



[21] 



9i-l 92 
Fi 



and 



From Eqs. ffTllTSllTTl) . it follows that 
and 



A [91+921 



Ai. = d{A2 - e^i) 



(22) 

(23) 

(24) 
(25) 
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where 



D 



l/[(l + 7V2)(2-i/)-2|/W/g2 



2/2(1 - 2mVQ2)(i + ^2) + 2(1 + i?)(l - y - 72?/2/4) 



d 



[1 + 7^y/2(l + 2m2y/g2)]^i _y_ ^2^2/4 
(l-2//2)(l + 722//2)-y2^2/g2 



r] = 7 



[l-|/-|/^(7V4 + m2/Q2)] 
-y/2)(l + 72y/2)-y2^2/g2 

1 - y/2 - y2m2/g2 



(26) 
(27) 
(28) 

1 +72y/2(l + 2m2y/g2)- ^^^^ 

Measurements of both A± and Ay yield the values of gi and 5^2 directly (given a parametrization 
of R and Fi), but in practice the majority of experiments in the past have only measured A\\. These 
experiments have been interpreted as measurements of gi, neglecting or correcting for g2. Although it 
turns out that g2 is small, there is no a priori reason for this, so it is safer to neglect some quantity for 
which we know an upper bound, e.g. we could utilize |y4_L| < 1, but a better approach is to utilize the 
(virtual) photon-nucleon asymmetries. The point is that there exists a restrictive bound on A2 (for the 
history of this see [HI [15] ) : 

1^2! < ^/R{l + Al)/2. (30) 



Replacing g2 in terms of gi and A2 yields 



All 



^ = (1+7- 



9^ 



+ (^ - 7) A2 



(31) 



where 77 — 7 is typically very small. Most experiments either use a parametrization for A2 or neglect 
the A2 terms in Eqs. ( !24l [3T|) altogether, leading to the approximations 



Ai^ — 
^ D 



;i+7^ 



9i_ 
Fi 



(32) 



to evaluate gi and then utilize the bound Eq. (l30l) to estimate the error involved. 

Strictly speaking, all of the above, from Eq. (JTj) onwards, applies only to a spin 1/2 target. For targets 
with higher spin, like the deuteron, the hadronic tensor contains many more structure functions. But 
in all the practical cases of interest the target nucleus may be considered as a weakly bound system of 
essentially independent nucleons, with a small correction for binding effects, and the nuclear asymmetries 
are then expressed in terms of the asymmetries on protons and neutrons. For a detailed analysis of this 
issue, and formulae relating measured asymmetries on deuterium and on '^He to the nucleon structure 
functions, see Section 2.1.4 of [12]. For the more general theory of DIS on spin 1 targets, see [16] and 

m 

Note that in the above we have kept terms of order M^/Q"^ and smaller. They are sometimes nec- 
essary in order to extract the correct experimental values of the structure functions from the measured 
asymmetries. However, very often, the QCD analysis of the structure functions is carried out at leading 
twist only, ignoring higher twist terms, i.e. of order M^/Q'^. This is clearly inconsistent in those cases, 
for example the Jefferson Lab experiments, where the above terms are important. 



1.2 Experiments 

When Bjorken proposed his famous sum rule [5] 40 years ago, he considered it "useless" since the 
experimental technology needed to test it appeared far from reach. It required substantial advances 
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in both polarized target and polarized beam technology before the first double-polarization lepton 
scattering experiment could begin [TS]. The first round of experiments at SLAC [TSl [121 [201 [21] using 
polarized electrons impinging on a polarized proton target seemed to confirm the expectations of the 
naive quark parton model (namely, that most of the proton spin was carried by its quark helicities). 
However, the range in x covered was rather limited, and the data had large errors and were taken at 
fairly low Q^. Nevertheless, these experiments started the exploration of the nucleon spin structure, 
not only in the DIS region, but also in the region of the nucleon resonances [20] . 

The EMC collaboration [2] used the (naturally polarized) muon beam at CERN, with much higher 
beam energy, together with a large polarized proton target, to extend the covered range down to 
significantly lower x ~ 0.01. Together with the earlier SLAC data, their result seemed to indicate that 
little to none of the proton spin was due to the helicities of the quarks. This violated the Ellis- Jaffe 
sum rule [3] within the simple quark-parton picture assuming that SU(3) is a good symmetry and the 
strange (sea) quark contribution to the nucleon spin could be ignored. The final data from the EMC 
experiment are published in [22] . 




Figure 2: Typical solid state polarized proton or deuteron target for electron scattering 
experiments. The cell containing the frozen ammonia (^^NHa or ^^NDs) is at the center of 
a Helmholtz-type magnet generating a homogeneous field of about 5 T. A ^He evaporation 
refrigerator (a liquid helium bath in a low-pressure environment) cools the target material 
to about 1 K. 140 GHz microwaves irradiate the target material to dynamically polarize the 
hydrogen nuclei. The polarization is measured by a resonant NMR circuit (the obtained 
signal vs. frequency is sketched in the top right). Polarized targets for muon beams are 
typically much longer and can be cooled to lower temperatures. 
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This puzzling result spurred the rapid development of several new experiments that had the goal to 
verify the data on the proton with greater precision and, most importantly at the time, test the Bjorken 
sum rule [S] by probing the spin structure of the neutron, as well. This last goal required the use of 
targets containing polarized neutrons, for which two very different technologies were developed. 

One possibility is to use polarized deuterons as a target composed of equal amounts of polarized 
protons and neutrons, from which neutron information can be extracted by comparison with pure proton 
targets (taking nuclear binding effects like the deuteron D-state into account). This method follows 
closely the example of polarized proton target technology, which typically involves chemical compounds 
(like alcohols or ammonia) rich in hydrogen (either for protons or = D for deuterons) which 
are seeded with paramagnetic centers (with unpaired electrons), cooled to cryogenic temperatures (1 
K or lower) and immersed in a strong magnetic field (usually several T). The unpaired electrons are 
polarized to nearly 100% in the field, and this polarization is transferred to the hydrogen nuclei via 
the process of Dynamic Nuclear Polarization (DNP), by employing microwave radiation of the proper 
frequency to induce coupled electron- nucleon spin transitions. This technique typically achieves a 
proton polarization of 80-90%, and a deuteron polarization of 30-50%, both of which can be measured 
using a Nuclear Magnetic Resonance (NMR) system. Figure [2] shows an example of such a target. In 
the case of deuteron targets, one can also use ^LiD compounds, where both the Lithium nucleus and the 
deuteron are polarized to increase the fraction of events from polarized nucleons (the so-called dilution 
factor) . 
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Figure 3: Typical gaseous ^He target for electron scattering experiments. The glass cell 
containing pressurized (~ 10 atmospheres) ^He gas is at the center of a Helmholtz magnet 
generating a homogeneous field of 25 gauss. The cell is connected to a second, heated 
volume where the ^He is polarized by spin-exchange collisions with alkali atoms, the latter 
being pumped by high-intensity lasers. The polarization is measured either via a NMR 
circuit (consisting of RF drive coils and NMR signal pickup coils) or using the EPR method. 



A complementary approach uses polarized "^He nuclei to gain information on the neutron, since to 
first order the "^He nucleus can be considered as a bound state of a proton-proton pair (with their spins 
anti-aligned) and an unpaired neutron carrying almost all of the nuclear spin. Again, nuclear corrections 
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to this picture are important, but the advantage is that one does not have to subtract a large contribution 
from polarized protons to extract the neutron results. The most widely used version of this target uses a 
glass vessel containing pressurized ^He gas mixed with a trace amount of an alkali atom vapor (typically 
Rubidium) that can be electron-spin polarized via optical pumping with circularly polarized laser light. 
The ^He nuclei are polarized through spin exchange collisions with the alkali atoms. A typical example 
is shown in Fig. [3l Polarizations of up to nearly 50% have been achieved in these targets, which can 
withstand much larger beam currents than the solid-state DNP targets. As an example, the polarized 
^He target used in JLab Hall A [23], with a 15 /xA electron beam on an 40 cm long, 10 atm ^He target, 
achieved a luminosity of 10'^^ cm^^s"^. Recently, a new development of a hybrid technique [21], using a 
mixture of Rubidium and Potassium atoms, has significantly increased the spin-exchange efficiency and, 
in turn, increased the target polarization up to 60% for the high density target. Further improvement to 
over 70% polarization has been achieved for this target by using recently available narrow-width lasers. 
The high luminosity and high polarization allow measurements of the spin structure functions with 
very high statistical precision. The target polarization is measured with both NMR and EPR (electron 
para-magnetic resonance) methods to a precision of 3%. Since the holding field is low, it is relatively 
easy to point the polarization in any direction in the horizontal plane using two pairs of Helmholtz 
coils. At Jefferson Lab, the addition of a new pair of vertical coils recently allowed the target also be 
polarized along the vertical direction, enabling the study of transverse spin. 

The successor collaboration to EMC at CERN, called SMC (Spin Muon Collaboration), used large 
dynamically polarized cryogenic deuteron [25] and proton [2B] targets to extract information on the 
neutron and to improve on the statistics and the kinematic reach of the EMC result. They also 
pioneered the use of semi-inclusive data, where a leading hadron is identified in coincidence with the 
scattered lepton, to get more information on the contribution of various quark fiavors to the nucleon 
spin [2Z]- The complete data set collected by the SMC resulted in precise inclusive results both at the 
highest momentum transfer [2S] and at the lowest quark momentum fraction x accessible to fixed 
target experiments [2U] . 

The E142 collaboration at SLAC was the first to use a '^He gas target with high luminosity to 
directly access the neutron spin structure functions and [5UI [5T| . Together with the EMC and 
SMC experiments, the results showed that the Bjorken sum rule including perturbative QCD (pQCD) 
corrections appeared to be valid. E142 was followed by a series of additional experiments at SLAC that 
used all three nuclear targets, proton, deuteron and ^He, to accumulate a highly precise data set on spin 
structure functions in the deep inelastic region. Instrumental for achieving ever higher precision was a 
significant improvement in the polarization (to over 80%) and intensity of available electron beams from 
strained GaAs cathodes irradiated with circularly polarized laser light. By using several electron beam 
energies and a set of up to 3 spectrometers, the E143 [321 ESI [M] and the E155 [35l [36] collaborations 
collected data on the proton and the deuteron over a wide range of momentum transfer at several 
values of x, which were used to study scaling violations for polarized structure functions. The E154 
collaboration added more neutron data at similar kinematics, using a polarized ^He target. The 

E143 collaboration also published the first precision results at lower and in the nucleon resonance 
region {W < 2GeV) [39]. The spin structure functions were measured with high precision by 

rotating the target polarization for all 4 experiments from a longitudinal to a perpendicular orientation 
to the beam [Si iQl E] . 

The most innovative approach to measuring DIS structure functions came from the HERMES col- 
laboration (see Fig. H]) which used positrons or electrons circulating in one of the HERA rings at DESY 
together with internal low-density gas targets fed directly from atomic beam sources [121 SS] . The tar- 
get atoms are polarized using hyperfine transitions induced by radio frequency fields and Stern-Gerlach 
type separation with magnetic sextupoles. The atomic beam is injected into a thin, windowless tube 
through which the beam circulates. This method yields a pure polarized target without any dilution 
from unpolarized materials. The polarization of the beam is accomplished by utilizing the Sokolov- 
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Figure 4: Layout of the HERMES experiment at HERA/Desy (Hamburg, Germany). The 
stored electron or positron beam traverses an open storage cell fed by an atomic beam source 
with polarized H, D or "^He. The scattered leptons (and leading hadrons) are detected in 
a large acceptance spectrometer consisting of wire chambers, scintillator hodoscopes, ring- 
imaging Cherenkov (RICH) and transition-radiation (TRD) detectors, and an electromag- 
netic calorimeter. 



Ternov effect (the spontaneous vertical polarization through spin-dependent synchrotron radiation of 
leptons in a storage ring). Spin rotators turn the polarization into the longitudinal direction at the 
target. The scattered electrons, as well as hadrons produced in coincidence, were detected by a large 
acceptance forward spectrometer. This setup allowed the HERMES collaboration not only to inde- 
pendently measure the inclusive spin structure functions gi and g2 (for final results see [S]), but also 
semi-inclusive structure functions for flavor-tagging [13] . In addition, they collected a large data set on 
related reactions of interest, from Deeply Virtual Compton Scattering (DVCS) [16] and transverse spin 
structure functions [17] to a first direct measurement of the gluon polarization [H] . 

After the shutdown of the HERA ring, there are now three laboratories left where experiments 
studying the spin structure of the nucleon continue: CERN, with the SMC-successor experiment "COM- 
PASS" ; BNL (on Long Island, NY) with the polarized proton-proton collision program at RHIC; and the 
Thomas Jefferson National Accelerator Facility ("Jefferson Lab" or JLab) in Newport News, Virginia, 
with an ongoing program of electron scattering in all 3 experimental halls. 

The COMPASS experiment (see Fig. [5]) uses the secondary (naturally polarized) muon beam at 
CERN together with large polarized deuteron and hydrogen targets to extend the kinematic reach and 
precision of SMC, and, as its main purpose, to extract information on the gluon polarization. This 
latter goal has been pursued by measuring both the production of hadron pairs with high transverse 
momentum [B] and by detecting charmed mesons in the final state (which are predominantly produced 
via photon-gluon fusion). Indirect information on the gluon contribution to the nucleon spin also comes 
from NLO analyses of inclusive DIS data, where the large kinematic lever arm offered by COMPASS 
makes an important contribution. First results on the deuteron have been published [l9l [501 EI] 
the COMPASS experimental program will continue in the foreseeable future. 

Another novel approach to studying the spin structure of the nucleon uses high-energy collisions of 
counter-circulating proton beams in the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Lab (BNL, Long Island, NY). Polarized protons are injected into a series of accelerators that finally fill 



11 




Figure 5: Layout of the COMPASS experiment at CERN (Geneva, Switzerland). 



both RHIC rings, where energies up to 100 GeV (250 GeV in the future) can be reached (see Fig. [6]). 
Siberian snakes rotate the proton spins to avoid depolarizing resonances, while spin rotators can select 
the desired spin direction at the interaction points. 

At present, there are two large experiments (PHENIX and STAR) that use polarized proton collisions 
to study the gluon helicity contribution AG to the nucleon spin. The observables used so far include 
meson production with high transverse momentum px as well as jet production, both probing the 
gluons through quark-gluon and gluon-gluon interactions in the initial state. First results from these 
experiments have been published ESI EH ES] and will be discussed in Section 12. 3[ By orienting the 
proton spins perpendicular to the beam direction, both experiments can also study reactions sensitive 
to transverse spins. 
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Figure 6: Layout of the RHIC accelerator complex at BNL (Long Island, NY). 
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Finally, for the last 10 years a large program using electron scattering to study the spin structure of 
nucleons has been underway at Jefferson Lab (JLab). This program has utilized the highest polarization 
electron beams (over 85%) with energies from 0.8 GeV to close to 6 GeV and all three species of polarized 
targets (p, d, and '^He) to study spin- dependent structure functions both in the DIS regime as well as 
in the nucleon resonance region (and even on nuclei). This program is ongoing in all three experimental 
halls and will be continued once the energy upgrade to 12 GeV of the JLab accelerator is completed in 
2014. In the following we give some of the experimental details for all three halls; the results achieved 
so far will be discussed in the relevant Sections l2.3f Hl 

The spin structure program in JLab's Hall A is focussed on the neutron, using a polarized '^He target 
(Fig. [3]) as an effective polarized neutron target. The target polarization direction can be oriented longi- 
tudinal or transverse to the beam direction. Measurements of polarized cross-sections and asymmetries 
in the two orthogonal directions allow a direct extraction of gi, g2, Ai and A2. A series of high precision 
experiments [561 EH EHl EHl [6OI [611 [621 [631 [Ml [65] measured gi and g2 in a wide range of kinematics, from 
very low Q'^ 0.01 GeV^) up to 5 GeV^ and from the elastic peak to the DIS region (W ^ 3 GeV). A 
pair of High Resolution Spectrometers (HRS) are used to detect the scattered electrons. The HRS have 
angular acceptances of ~ 6 msr and momentum acceptances of ~ 9%. Their angular range is 12.5° to 
160° and can reach as low as 6° with the addition of a septum magnet. The high luminosity of 10^^ 
cm~^s~^ allowed for precision measurements at numerous HRS momentum and angular settings to cover 
a wide swath in the (Q^,W)-plane. The electron detector package consists of vertical drift chambers (for 
momentum analysis and vertex reconstruction), scintillation counters (for data acquisition trigger), gas 
Cherenkov counters and lead-glass shower calorimeters (for particle identification). The HRS optical 
property and acceptance have been carefully studied. Absolute cross sections are measured to a level of 
2-3% precision. Asymmetries are measured to a level of 4-5% precision, mostly due to the uncertainties 
from the beam and target polarization measurements. The spin structure functions gi and (72 are ex- 
tracted using polarized cross section differences in which contributions from unpolarized materials, such 
as target windows and nitrogen, cancel. Corrections for the two protons in '^He are still needed since 
they are slightly polarized due to the D state ( 8%) and S' state (1.5%) of the ^He wave function [66] . 
Corrections for binding and Fermi motion are applied using state-of-the-art ^He calculations [671 168] . 
Uncertainties due to the nuclear corrections have been studied [67]. In the region of DIS and for the 
extraction of moments, the uncertainties are usually small, typically less than 5%. 

The EG1-EG4 series of experiments in JLab's Hall B (see Fig. [7]) has as its goal to map out 
the asymmetry Ai and the spin structure function gi of both nucleons over the largest, continuous 
kinematic range accessible. It makes use of the CEBAF Large Acceptance Spectrometer (CLAS) in 
Hall B that covers an angular range of about 6 degrees to over 140 degrees in polar angle and nearly 
27r in azimuth [69]. Because of the open geometry and the toroidal magnetic field (maintained by 6 
superconducting coils evenly distributed in azimuth), one can simultaneously detect scattered electrons 
over a wide kinematic range, as well as secondary produced hadrons (nucleons, pions and kaons) for 
semi-inclusive or exclusive channels. By combining runs with several different beam energies from 1 to 
6 GeV, a continuous coverage in from 0.015 to 5 GeV^ and in final state mass W, from the elastic 
peak {W = 0.94 GeV) to the DIS region {W ~ 3 GeV), has been achieved. Inclusive results from the 
EGl experiment have been pubhshed [701 \7T\ [72l [73] . 

The relatively low luminosity limit of this open detector (about 10^'* cm~^s~^) is rather well matched 
to the luminosity limits typical for solid state, dynamically polarized targets. EG1-EG4 took data on 
both hydrogen (^^NHs) and deuterium (^^NDs). So far, only targets polarized along the beam direction 
have been utilized (because of the difficulty to combine a large transverse magnetic field with the 
open geometry of CLAS), which necessitates (minor) corrections of the measured asymmetries for the 
unobserved contribution from A2 (see Eg. |3T1) . A fit to the world data on A2 and on unpolarized structure 
functions R and Fi [7^ [75] is used to extract the desired spin structure function information from the 
measured asymmetries. In addition to the structure function gi{x, Q"^), the CLAS data have also yielded 
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Figure 7: Schematic of the CEBAF Large Acceptance Spectrometer (CLAS) with the polar- 
ized target. Charged particle tracks are measured with three layers of drift chambers, while 
electrons can be identified using a set of Cherenkov counters, time-of-flight scintillators and 
a electromagnetic calorimeter (the latter also shown in end-on view). 



new results on resonance excitation and decay (via exclusive ir^, ir^ and vr" channels) [761177] . on deeply 
virtual Compton scattering [78j, and on single and double spin asymmetries in semi-inclusive hadron 
production [79] . 

The first experiment completed in Hall C used a standard DNP ammonia target (^^NHa and ^^NDa) 
and the existing high momentum spectrometer (HMS) for a detailed look at the resonance region at 
intermediate ~ 1.3 GeV^. This is the only experiment at JLab on the proton and the deuteron where 
both longitudinal and transverse double spin asymmetries were measured, allowing an unambiguous 
separation of the structure functions Ai and A2 or gi and (72 up to a final state missing mass of TV ~ 2 
GeV. First results have been published [80] . 

2 Spin-dependent Parton Density Functions 

2.1 The Simple Parton Model 

The parton model was invented [81] long before QCD to explain the fact that the unpolarized structure 
functions Fi^2 apparently scale, i.e., do not decrease with Q^, in contrast to the elastic form factors 
which decrease rapidly [1] . The partons were conceived of as effectively massless, point-like constituents 
of the nucleon, which interact electromagnetically like leptons, and it was argued that a nucleon, in 
a frame where it is moving very fast, could be viewed as a "beam" of coUinear partons, as shown in 
Fig. [HI The partons are characterized as having momentum p = x'P and covariant spin vector s. 

The interaction with the hard photon is then visualized as in Fig. [9l in which the lepton-quark 
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P.S 




p- = x'. P, s. 



Figure 8: The nucleon viewed as a coUinear beam of partons. 

scattering is treated analogously to elastic lepton-lepton scattering. Requiring the final quark to be on 
mass shell, i.e. (p + g)^ = 0, selects the value x' = x. 



lepton 



nucleon 
*- 




Figure 9: Parton model description of DIS. 

For unpolarized DIS one finds the scaling result expressed in terms of the number density q{x) of 
quarks and q{x) of antiquarks 

Fiix,Q') = lj2'''M^) + ^A^)] (33) 
j 

where the sum is over fiavors j, ej is the charge of the quark, and the Callan-Gross relation [82] yields 

F2{x) = 2xFi{x). (34) 
For longitudinally polarized DIS one obtains 

j 

with 

Aq{x) = — (36) 

where q±{x) are the number densities of quarks whose spin orientation is parallel or antiparallel to the 
longitudinal spin direction of the proton (see Fig. [TOl) . In terms of these, the usual (unpolarized) parton 
density is 

q{x) = q+{x) + q-{x). (37) 




Figure 10: Visualization of the longitudinally polarized parton density Aq{x). The upper 
arrows show the spin direction. 

Concerning the second spin structure function, there are many different, inconsistent results for 
g2{x) in the literature. In fact there is no unambiguous way to calculate g2 in the simple parton model. 
For an explanation of the problem see Section 3.4 of [12]. 
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The only reliable result is the Wandzura-Wilczek relation [9] 



92[X) 



X 



dx' 



(3J 



which was originally derived as an approximation in an operator product expansion approach, but 
which has recently been shown to be derivable directly in the simple parton model |83j . It is discussed 
further in Sections 12.61 and [31 

Let us rewrite the expression for gi, Eq. fl5^ . in terms of linear combinations of quark densities 
which have specific transformation properties under the group of flavor transformations SU{?))f'- 



Aqs = (An + An) - (Ad + Ad) 

Aqs = (An + An) + {Ad + Ad) - 2{As + As) 
AS = (An + An) + {Ad + Ad) + (As + As). 



(39) 

(40) 
(41) 



These transform respectively as the third component of an isotopic spin triplet, the eighth component 
of an SU{3)f octet, and a flavor singlet. Then 



^Aq3{x) + ^Aqs{x) + AE{x) 



(42) 



Taking the first moment of this yields, for protons. 



^1 = J 9i{x)dx = - 



3 1 

-as + -as + ao 



(43) 



where 



as 



dx Ag3(x) 



dx Aqs{x) 



ao 



'0 

AS 



dx AS(a;) 



(44) 



Now the hadronic tensor W'^'^ is given by the Fourier transform of the nucleon matrix elements of 
the commutator of electromagnetic currents Jn{x): 



W,,{q;P,S) = ^ f d\ e'^-%P,S\[J,{x),JM]\P.S). 
Ztt J 



(45) 



In hard processes, ~ is important, so we can use the Wilson Operator Product Expansion (OPE). 
This gives moments of gi^2 in terms of hadronic matrix elements of certain operators multiplied by 
perturbatively calculable Wilson coefficient functions. This is discussed in more detail in Section 13.11 
The in Eq. (jH]) turn out to be hadronic matrix elements of an octet of quark SU{3)f axial- vector 

The octet of currents 



currents Jg^ {j 



1, 8) and a flavor singlet axial current J^^ (see Eqs. 



is precisely the one that controls the weak decays of the neutron and of the octet hyperons, which implies 
that the values of 03 and ag are known from /3-decay measurements: 



ag = c/^ = 1.2670 ±0.0035 ag = 0.585 ± 0.025. (46) 
Hence a measurement of Fi, Eq. ( 143|1 . can be considered as giving the value of the flavor singlet Oq. 
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It should be noted that the connection with hadronic matrix elements of local operators is strictly 
only valid if the moments include the elastic contribution, which appears as a delta-function at a; = 1. 
However in deep inelastic scattering the elastic contribution is completely negligible, so the experi- 
mentally tabulated moments do not include an elastic term, as can be seen in the discussion of the 
extrapolation to x = 1 given in Sections 12.41 and 12. 5[ However, in dealing with inelastic scattering at 
low Q^, especially in the resonance region, it is important to distinguish between moments including the 
elastic contribution, for which the notation will remain F and those not including an elastic contribution 
which will be denoted F. 

Now we can put the startling result of the European Muon Collaboration [2] into context: Knowing 
the values of as and as, the EMC measurement implied 

a^*'^^ ^ 0. (47) 

But in the naive parton model 

ao = AS = as + 3(As + As) (48) 

where AS is given by Eq. (HTj) . 

In 1974 Ellis and Jaffe [3] had suggested that one could ignore the contribution from the strange 
quarks, i.e., from As -|- As, implying that 

ao ~ as ^ 0.59. (49) 

Thus the EMC result Eq. (147!) is in gross contradiction with Ellis- Jaffe. It was this contradiction which 
at first aroused interest in the EMC result, but it was soon realized that their result had far more 
serious consequences. 

Consider the physical significance of AS(a;). Since q±{x) count the number of quarks of momentum 
fraction x with spin component ±| along the direction of motion of the proton (say the z-direction), 
the total contribution to Jz coming from the spin of a given flavor quark is 

Jo L v^y V ^ y j 

(50) 

^ Jo 
It follows that 

ao = 2{Sf"'''^') (51) 

where (^Sf^"''^'^^) is the contribution to from the spin of all quarks and antiquarks. 

Naively, in a non-relativistic constituent model one would have expected all of the proton spin to 
be carried by the spin of its quarks. In a more realistic relativistic model one expects 2(S'f"'^''^*) ~ 0.6, 
which, as first noted by Sehgal [84], is actually close to the value obtained by neglecting the strange 
contribution in Eq. (HSI) . but quite far from the EMC value Eq. (H7I) for a^. 

It was this discrepancy between the contribution of the quark spins to the angular momentum of the 
proton, as measured in DIS and as computed in both non-relativistic and relativistic constituent models 
of the proton, that was termed a "spin crisis in the parton model" [1] . It should be noted, however, that 
Thomas and co-workers have developed a constituent model, a more sophisticated version of the Cloudy 
Bag model, in which the nucleon is visualized as made up of valence quarks and a pion cloud, and the 
wave function includes the effect of single gluon exchange between the valence quarks. The pion cloud 
and the strength of the one gluon exchange were treated phenomenologically, determined from fits to 
various different static properties of the hadron octet and decuplet, and it was not clear if , or to what 
extent, double counting was taking place. However, lattice studies of the mass-splitting between the A 
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and the nucleoli [85] suggest that the pion cloud contributes very little to the splitting and thus allows 
a fairly reliable estimate of the value of the phenomenologically determined gluon coupling strength. 
It then turns out that the contribution of the quark spins to the angular momentum of the proton, in 
its rest frame, is AS ^ 0.35, which is not too far from most recent experimental results. Thus in this 
approach there is no significant "spin crisis". For a review of this work see |86j . 



2.2 The Parton Model m QCD 

The parton model is a heuristic picture, much like the impulse approximation in nuclear physics, and 
predates QCD. Once QCD is accepted as the theory of strong interactions, with quark and gluon fields 
as the fundamental fields, there will be interaction dependent modifications of the simple parton model 
formulae for DIS. The theory is invariant under color gauge transformations, but the physical content of 
individual Feynman diagrams does depend upon the gauge, and it turns out that a parton-like picture 
emerges in the light-cone gauge Aj^ = where is the gluon vector potential. The description of 
nucleon structure becomes much more complicated, involving a set of twelve functions, and the parton 
model number densities q{x), Aq{x) (and the analogue for transversely spinning nucleons Aj'g(x), 
which is not relevant for the longitudinal spin structure) are only the principal, so called "leading 
twist" members of this set. In this review, we shall only deal with q{x) and Aq{x). The reaction is 
visualized as in Fig. [TT], where the top "blob" involves a hard interaction (the photon is highly virtual) 
and the bottom "blob" involves non-perturbative soft interactions. 




Figure 11: QCD generalization of parton model. 



The main impact of the QCD interactions will be twofold: 

1) to introduce a mild, calculable logarithmic dependence in the parton densities 

2) to generate a contribution to gi arising from the polarization of the gluons in the nucleon 

We shall not go into the technical details, but simply indicate the physical source of these effects. 

1) QCD corrections and evolution. 
In Fig. [12] we show the Born term for the interaction of the virtual photon with a quark (the hard blob), 
and the simplest correction terms, a vertex correction and a diagram where a gluon is radiated from 
the active quark before it interacts with the photon. 

Unfortunately these correction terms are infinite. The infinity is caused by coUinear divergences 
which occur because of the masslessness of the quarks and which are removed by a process known 
as factorization. In this process the reaction is factorized (separated) into a hard and soft part and the 
infinity is absorbed into the soft part, which in any case cannot be calculated and has to be parametrized 
and studied experimentally. 

The point at which this separation is made is referred to as the factorization scale fi^. Schematically, 
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(b) 



Figure 12: Example of QCD correction terms (b), to the Born approximation (a). 



one finds terms of the form In ^ which one sphts as follows 



Oig In — = as In + In 




(52) 



and one then absorbs the first term on the right hand side into the hard part and the second into 
the soft part. yU.^ is an arbitrary number, like the renormalization scale, and, in an exact calculation, 
physical results cannot depend on it. However it does mean that what we call the parton density has 
an extra label /i^ specifying our choice. Moreover, since we never calculate to all orders in perturbation 
theory, it can make a difference what value we choose. It turns out that an optimal choice is /i^ = Q^, 
so the parton densities now depend on both x and i.e. we have q{x, Q"^) and Aq{x, Q^), and perfect 
Bjorken scaling is broken. But the variation with is gentle (logarithmic), and can be calculated via 
what are called the evolution equations which will be discussed later. 

It turns out to be crucial in handling these divergences to use the technique of dimensional reg- 
ularization, which is straightforward in the unpolarized case, but which runs into ambiguities in the 
polarized case. As a consequence there are several different factorization schemes in use and it is crucial, 
when presenting results on the parton densities to specify which scheme is being utilized. 

At present there are three schemes in use: 

i) The Vogelsang, Mertig, van Neerven scheme [871 EH], MS — MNV (usually just abbreviated as 
MS). In this scheme as and ag are independent of Q^. 

ii) The AB scheme of Ball, Forte and Ridolfi [89], which, in addition, has the first moment 




independent of Q^. 

iii) The JET scheme of Carlitz, Collins and Mueller [90], Anselmino, Efremov and Leader [12] and 
Teryaev and Miiller [91], and which is identical to the Chiral Invariant scheme of Cheng [92]. In this 
scheme 03 and as are independent of as is AS, but it can be argued that the JET scheme is superior 
to the others in that all hard effects are included in H^i, (see Fig. [TTl) . 

Of course, if one could work to all orders in perturbation theory it would make no difference which 
scheme one used, but given that we work to leading order (LO), next to leading order (NLO), and in 
some cases to NNLO, the choice of scheme can be of importance. 





(53) 
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For the polarized densities the evolution equations are 
d A / ctsiQ^) dy 



dlnQ^ 



Aq{x,Q') = ^4P- I ^{APMy)Aq{y,Q') 

AP,Gix/y)AG{y,Q')} (54) 



dlnQ 

where AG{x) is analogous to Aq{x 



AG(x,Q2) = ^ r ^{APa,ix/y)Aqiy,Q') 

+ APGG{x/y)AG{y,Q')}, (55) 



AG(x) =G+(x) -G_(x). (56) 
The AP are the polarized splitting functions and are calculated perturbatively 

AP{x) = AP(°) (x) + — AP(^) (x) (57) 

27r 

where the superscripts (0) and (1) refer to LO and NLO contributions. For details about these the 
reader is referred to [87]. A pedagogical introduction to scheme dependence can be found in Section 
11.8 of [93]. For a more complete set of relations between the densities and Wilson coefficients in the 
various schemes, see [M] . 

Note that in LO flavor combinations like qj — q^ (e.g., u{x) — d{x)) and valence combinations like 
Qf ~Qf (s-g-, u{x) — u{x)) are non-singlet and evolve in the same way, without the AG term in Eq. fl5^ . 
(There is no splitting in LO from a g to a g, nor from, say, a m to a d.) However, in NLO flavor 
non-singlets like u{x) —d{x) and charge-conjugation non-singlets like u{x) —u{x) evolve differently. The 
origin of this difference can be seen in Figs. [I3] and [TH Fig. [13] shows an NLO amplitude for a quark to 
split into a q . 




Figure 13: NLO amplitude for the q q transition. 

Fig. [T^ shows two possible contributions to APgg from taking the modulus squared of this amplitude. 
In (a) the contribution is pure flavor singlet and involves only gluon exchange, whereas in (b) the 
contribution is non-singlet. However, if we try to do something similar for a flavor changing splitting 
function e.g. APdu we find that we cannot construct the non-singlet diagram. 

The expression for gi{x,Q'^) now becomes 

9i{x,Q') = ^ E el{Aq{x,Q') + Aq{x,Q') 

flavors 

asm 



+ 



27r 



"^{ACgix/y) [Ag(y, Q') + Ag(y, Q')] 

J X y 



ACG{x/y)AG{y,Q')}] (5J 
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(a) (b) 
Figure 14: NLO contributions to the splitting function for a. q —>■ q transition. 



where ACq and ACq are Wilson coefficients evaluated from the hard part calculated beyond the Born 
approximation. Note that very often the evolution equations are written using the convolution notation, 
for example, 

AC,®Ag= / ^ACqix/y)Aqiy). (59) 

J X y 

2) The gluon contribution to gi. 
Fig. [T5l demonstrates a NLO gluon-initiated contribution to DIS. 




Figure 15: QCD diagram leading to the anomalous gluon contribution. 




k, A 



k, A 



Figure 16: Feynman diagram responsible for the anomaly. 



In the Bjorken limit, for the longitudinal polarized case, this turns out to involve the gluonic version 
of the famous Adler [95] , and Bell and Jackiw [96] anomalous triangle diagram shown in Fig. [161 The 
net result is that there is an anomalous gluonic contribution to the flavor singlet ao [901 IHZl IMl 199] 



a: 



gluons 




= -3 
= -3 



asm 
asm 

2ti 



dx AG{x, 



AG{Q^). 



(60) 
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Note that the factor 3 corresponds to the number of light flavors i.e. u, d, s. Heavy flavors do not 
contribute. 

So, there exists potentially a gluonic contribution to the first moment of gi. 

rf """^(g^) = _i ^^ii^ AG(g2) i (ei) 

3 27r 

This result is of fundamental importance. It implies that the simple parton model formula Eq. (jHj) for 
ao (and hence for F^) in terms of the Aqf could be considered as incomplete. Instead, 

ao = AS - 3 AG. (62) 

However, it should be noted that there is a subtlety concerning Eq. fl62l) . It turns out that the 
result actually depends on the factorization scheme utilized. Eq. ( l62i) is correct in the AB and JET 
schemes, but the gluon contribution to oq is zero in the MS scheme. Now the spin crisis emerged from 
comparing the size of AS with what would be expected if the quark spins dominated the nucleon's 
angular momentum. But in NLO, in the MS scheme, AS varies with Q^, so one could argue that it 
cannot be directly interpreted as the spin contribution of the quarks. It is the invariant AB or JET 
value of AS which should be interpreted directly as a spin. 

The fundamental conclusion is that the small measured value of Oq does not necessarily imply that 
the physically meaningful, invariant AS is small. This discovery was hailed as a resolution of the spin 
crisis. We shall see later that this seems to be a false hope. 



2.3 Experimental Determination of Polarized Parton Densities 
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Figure 17: World data on the polarized structure function xgi{x) for the proton, deuteron 
and neutron in the DIS region {W > 2 GeV), taken by different experiments, at several 
different values of Q^, as compiled by the Particle Data Group [100] (left panel). The proton 
data are also shown versus Q^, for several bins in x |101] (right panel; the dashed lines 
just connect the data to guide the eye). Note that much of the data at very small x is for 
Q2 < 1 GeV^. 
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From the discussion in the previous two sections, it is clear that the bulk of information on polarized 
parton densities in the nucleon comes from charged lepton scattering experiments in the deep inelastic 
(DIS) region, with final state masses larger than W = 2 GeV and momentum transfers in excess of 

= 1 GeV^. As described in Section ll.2[ a vast amount of data on the inclusive spin structure 
function gi{x, Q^) has been accumulated by lepton scattering experiments at SLAG, GERN, DESY and 
Jefferson Lab. The Particle Data Group's compilation |100] of these data plotted vs. x are shown in 
Fig. [171 on the left side, while the right side displays the Q^-dependence of gf{x, Q"^) for several fixed x 
bins. The observed mild violation of scaling is sensitive to the polarized gluon density (see Sections 12.21 
and l2.5p . In Section [231 we highlight in particular the new data taken at high values of x, where valence 
quarks dominate the measured spin structure functions. In that region, valence parton densities can be 
uniquely determined from DIS data on the proton and the neutron. 




Figure 18: World data on polarized parton densities Aq[x) extracted from semi-inclusive 
DIS data, as compiled by the Particle Data Group |1UU] . 



Inclusive data alone, however, are not able to distinguish the contributions from quark and antiquark 
densities. One approach to gather additional information has been to use semi-inclusive lepton scattering 
(SIDIS), where in addition to the scattered lepton one detects a leading hadron (typically a pion or 
kaon) in the final state. Struck quarks of different flavors have varying probabilities to produce those 
final state hadrons, as expressed through the various fragmentation functions. These are principally 
extracted from e^e~ collision experiments, and still have considerable uncertainties (in particular for 
kaons). By comparing the asymmetries of semi-inclusive 7r"'",7r~ and K~^,K~ production from protons 
and neutrons, one can try to separate out the contributions from all quark flavors. This approach 
was pioneered by the SMG collaboration [27] and the most detailed data set stems from the HERMES 
collaboration [351 I102j . as shown in Fig. [181 The GOMPASS collaboration has also measured semi- 
inclusive channels [51] . 
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While the virtual photon in deep inelastic lepton scattering couples directly to the quarks (via their 
electric charges), gluons can only be probed indirectly in this reaction. Their main influence is in the 
gentle (logarithmic) evolution of the quark densities with (see r.h.s. of Fig. [T7I) . as outlined in 
Section I2.2[ The resulting violation of Bj or ken-scaling can be exploited to extract the contribution 
of gluons to the nucleon momentum and spin. This has been spectacularly successful in the case of 
unpolarized DIS, where (thanks to the existence of the lepton- nucleon collider HERA) the range of 

that has been accessed experimentally exceeds 5 orders of magnitude for some (fixed) values of 
X. The situation is much more difficult in the case of polarized structure functions; so far, only fixed 
target experiments have measured double spin observables in lepton-nucleon scattering. The highest 
available are of order 100 GeV^ (using the muon beam at CERN). Recent data from COMPASS [19] have 
dramatically improved the precision of g\ at relatively high for the deuteron which is an isoscalar 
target and therefore more sensitive to gluon contributions. For the maximum lever-arm in Q^, it is 
important to push for high precision at both extremes of the range. In this context, new deuteron 
results from Jefferson Lab [72] play an important role. 
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Figure 19: Results (some preliminary) on the gluon polarization AG/G extracted (at leading 
order) from photon-gluon fusion processes. (Data compiled by the COMPASS collabora- 
tion [T03] ). 



For a direct approach, it is highly desirable to measure observables that are more closely linked 
to the polarized gluon density in the nucleon, e.g. virtual or real photon-gluon fusion (see Fig. [15!) 
in semi-inclusive lepton scattering or photon absorption, with two quarks in the final state stemming 
from the splitting of a gluon. In order to distinguish this process from the more copious direct quark 
diagrams, one has to select either events in which the final state hadrons or jets carry large transverse 
momentum "px or where the final state quark has a heavy flavor (e.g., charm) unlikely to stem from the 
nucleon itself. This approach has been used by HERMES and COMPASS and was to be employed by 
the (unfortunately terminated) SLAC experiment E161. 

The existing results (including preliminary data) are shown in Fig. [191 The new Hermes data 
point (inverted triangle) is extracted from high-p^- virtual photon production of a single hadron, while 
the COMPASS data points come from high-p-p hadron pairs or jets and from open charm production 
data [jy^ and D*; indicated by the star symbol in Fig. [T9ll . These data put some restrictions on the 
magnitude of the polarized gluon density AG(x) in the moderate x-region to which they are sensitive 
but still have rather large statistical errors. Also, at present no NLO analyses of these data exist, and 
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backgrounds from "ordinary" SIDIS processes as well as "resolved photon" contributions have to be 
taken into account. 

An alternate route has been employed in the polarized proton collision program at BNL, using the 
Relativistic Heavy Ion Collider (RHIC). One selects final state signatures like hadrons, jets or even 
direct photons with large transverse momentum that indicate an underlying hard interaction between 
two constituents from the two colliding protons. The cross section for the production of this final state 
can be expressed as a convolution of the parton densities of the two initial protons with the hard parton- 
parton scattering cross section and the final state fragmentation function (see, e.g., |104] ). There is 
good agreement between the NLO pQCD predictions and data for the unpolarized cross section |54j . 
suggesting that pQCD can be applied reliably to such reactions. 
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Figure 20: Double spin asymmetries in high "px 7r° production (PHENIX [105] , left panel) 
and jet-production (STAR preliminary [106] , right panel) from proton proton collisions at 
200 GeV center-of-mass energy at mid-rapidity. The various curves shown are the NLO 
pQCD predictions from several versions of the polarized gluon density parametrizations 
from Vogelsang et al. |107] . 



For the double spin asymmetry, one divides the cross section difference for opposite proton spins by 
the sum. The result depends on products of polarized parton densities; depending on the kinematics 
one is mostly sensitive to quark-quark, quark-gluon or gluon-gluon subprocesses. While this quadratic 
dependence does not allow a "naive" direct extraction of polarized parton distributions (not even in 
LO), the observed asymmetry can be compared to predictions from parametrizations of quark and gluon 
densities. Figure 1211] shows the preliminary results from PHENIX (vr^ production, |105j ) and STAR (jet 
production, [106] ) for the double spin asymmetry Aj^i^ versus transverse momentum p^, compared with 
several model calculations based on the GRSV parametrization |108] of polarized parton densities. One 
clearly can see a preference for the curves labeled "AG = 0" while some of the more extreme possibilities 
(AG = ±G) can be excluded. However, depending on the particular overall shape of the polarized gluon 
distribution with x, non-zero values of the overall first moment of AG are not excluded. 

The most comprehensive analysis of all DIS, SIDIS and pp data so far has been conducted by de 
Florian, Sassot, Stratmann and Vogelsang |104] . The authors find good agreement between all data 
sets, which together yield a preference for the integral over AG(a;) from x = 0.05 to a; = 0.2 to be close 
to zero. For a definite answer on the question "what is the contribution AG to the nucleon spin" we 
need additional data that can better constrain the higher and lower x region, to pin down both the 
shape and overall magnitude of AG(x). These data will be forthcoming from future experiments at 
JLab (precision measurements of DIS on the deuteron, both at 6 GeV and at 12 GeV beam energy) and 
the continuation of the COMPASS experiment at CERN, as well as higher precision data for different 
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center-of-mass energies and different final state channels from RHIC. 

In Section 12.5^ we survey the status of the polarized parton densities extracted from all the data 
described above. 

2.4 Measurements of Spin Structure in the Valence Region 

The properties of nucleon structure functions in the high-x region are of special interest, because this is 
where the valence quark contributions are expected to dominate. With minimal sea quark and explicit 
gluon contributions, it is a clean region to test our understanding of nucleon structure. Relativistic 
constituent quark models |1U91 [SB] should be applicable in this region and perturbative QCD [IIU] can 
be used to make predictions in the large x limit. 

To first approximation, the constituent quarks in the nucleon are described by SU(6) wave functions 
with orbital angular momentum of zero. SU(6) symmetry leads to the following predictions jlllj : 

d{x)/u{x) = 1/2; Au{x)/u{x) = 2/3; Ad{x)/d{x) = -1/3; A{{x) = 5/9; and = 0. (63) 

Relativistic Constituent Quark Models (RCQM) with broken SU(6) symmetry, e.g., the hyperfine 
interaction model |1U9] . lead to a dominance of 'quark-di quark' configurations with the spectator- 
diquark spin S" = at high x. This implies that clS X ^ 1 '. 

d/u 0; Au/u 1; Ad/d -1/3; Al 1; and 1. (64) 

In these RCQM models, relativistic effects lead to a non-zero quark orbital angular momentum and 
reduce the valence quark contributions to the nucleon spin from 1 to 0.6 — 0.75. 

Another approach is leading-order pQCD |110j . which assumes the quark orbital angular momentum 
to be negligible and leads to hadron helicity conservation. For a nucleon with helicity +1/2 

g±(x) ^ (l-x)2'^-i+(^=Fi) (g5) 

as a: 1, where n = 2 is the number of spectators. The exponent for q_^. is therefore 3, while the 
strength of g_ is suppressed by an extra factor of (1 — x)^ at large x. This, plus SU(6), leads to: 

d/u 1/5; Au/u 1; Ad/d ^1; ^ 1; and A^ 1. (66) 

in agreement with [112j . Eq. fl65l) cannot hold at all Q^, since evolution causes the power of (1 — x) to 
grow like loglogQ"^ at large Q^, but Eq. fl66|) should hold at all Q/^. 

Not only are the limiting values at x = 1 important, but also the behavior as x approaches 1, which 
is sensitive to the dynamics in the valence quark region. A new approach |113] in pQCD, including 
quark orbital angular momentum (and therefore not requiring hadron helicity conservation), shows a 
different behavior at large x region while keeping the same limiting values at x = 1. In particular, the 
approach of Ad / d towards unity only sets in at a significant larger x than predicted by hadron helicity 
conservation (see below). 

Experimentally, it is difficult to access the high-x region because most high-energy experiments so 
far lacked the necessary luminosity and fine-enough resolution in x. This situation has dramatically 
changed with the advent of the high-current, moderate energy continuous electron beam at Jefferson 
Lab. 

JLab Hall A experiment E99-117 [HUES] measured the neutron asymmetry A" with high precision 
in the x region from 0.33 to 0.61 (Q^ from 2.7 to 4.8 GeV^). Asymmetries from inclusive scattering 
of a highly polarized 5.7 GeV electron beam on a high pressure (> 10 atm) (both longitudinally 
and transversely) polarized ^He target were measured. Parallel and perpendicular asymmetries were 
extracted. After taking into account the beam and target polarizations and the dilution factor, they 
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Figure 21: (left-panel) and A? (right-panel) results from JLab Hall A E99-117 [6ll[62] and 
CLAS EG lb [72] experiments (filled circles), compared with the world data and theoretical 
predictions (see text for details). 



were combined to form . Using a recent model [68], nuclear corrections were applied to extract A". 
The results on are shown in the left panel of Fig. [2T1 The experiment greatly improved the precision 
of data in the high-x region, providing the first evidence that A^ becomes positive at large x, a clear 
evidence for SU(6) symmetry breaking. The results are in good agreement with the LSS 2001 pQCD 
fit to previous world data |114] (solid curve) and the statistical model |115j (long-dashed curve). The 
trend of the data is consistent with the RCQM |109] predictions (the shaded band). The data disagree 
with the predictions from the leading-order pQCD models [110] (short-dashed and dash-dotted curves). 
These data provide crucial input for the global fits to the world data to extract the polarized parton 
densities and the extractions of higher-twist effects. 

New results for A\ and Af from the Hall B EGlb experiment [72] are also available. The data cover 
the range of 1.4 to 4.5 GeV^ for x from 0.2 to 0.6 with an invariant mass larger than 2 GeV. The 
results on A^ are shown in the right panel of Fig. [2IJ The precision of the data improved significantly 
over that of the existing world data. Similar data also exist on the deuteron and once again exhibit a 
trend to exceed the asymmetries predicted by SU(6) (Eq. [63]) at large x. In addition to the prediction 
by the RCQM model [lOQj (again indicated by the shaded band), several curves based on different 
scenarios of SU(6) symmetry breaking as presented in the paper by Close and Melnitchouk [116j are 
also shown in the right panel of Fig. [2TJ These scenarios are based on duality arguments and include 
final state helicity-1/2 dominance (dashed curve), spin-1/2 dominance (dotted) and symmetric wave 
function suppression (dash-dotted), which is closest to the data. 

The polarized quark distribution functions Au/u and Ad/d in the high-a; region were first extracted, 
in the leading-order approximation, from the Hall A neutron data and the world proton data (left panel 
of Fig. [221) [611 |62] . A recent leading-order extraction used the CLAS EGlb proton and deuteron 
data [72] as well. The results are shown in the right panel of Fig. [22], along with predictions from 
leading-order pQCD [110] (dashed curves) and a pQCD fit (solid line) including quark orbital angular 
momentum contributions [113] . The results of Ad/d are in significant disagreement with the predictions 
from the leading-order pQCD model assuming hadron helicity conservation. Data agree better with the 
fit including quark-orbital angular momentum contributions, suggesting that the quark orbital angular 
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Figure 22: Au/u (upper side of both panels) and Ad/d (lower side of both panels) results 
from JLab Hall A E99-117 [HH [62] (left panel) and from CLAS EGlb [72] (right panel), 
compared with the world data and theoretical predictions (see text for details). 



momentum may play an important role in this kinematic region. 
2.5 Status of Polarized Parton Densities 

From measurements of gi{x,Q^) at various values of x and Q^, one can, in principle, via Eq. (158|) . 
extract the quark and gluon polarized parton densities. We shall first consider general issues related to 
the determination of the parton densities and then discuss what is known about the densities themselves. 

1) General properties of the parton densities 

• In DIS, one can only obtain information on the combinations Aq{x, Q^)+Ag(x, Q^), i.e., there is no 
information at all about the Am and Ad sea quarks, though As + As is determined. Nonetheless, 
in presenting results, some authors use conventions such as Am = Ad = As, or variants thereof, 
in order to show results for the valence quarks. It is important to bear in mind that these results 
are not physical and are entirely convention dependent. 

• To determine Am and Ad one needs to study semi-inclusive deep inelastic scattering or semi- 
inclusive hadron-hadron reactions, where, however, one also needs information about fragmenta- 
tion functions. 

• Because the gluon has zero electric charge it does not couple directly to the photon, so it con- 
tributes in Eq. flSSl) only in NLO. However, as shown in Eq. ([SID, ^l^o plays a role in the 
evolution of the parton densities, but since the evolution is only logarithmic and the lever arm 
in Q'^ of present day experiments is relatively small {much smaller than in the unpolarized case), 
there is considerable uncertainty about the form and magnitude of AG{x, Q"^). 

• The parametrizations assumed for the Aq{x) at some value of follow general intuition based 
on the parton model. In particular, if the q+-{x) are regarded as positive number densities then, 
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via Eqs. fl36| [37j) one imposes the positivity condition 

\Aq{x)\ < q{x). (67) 

It can be shown that if positivity holds at some Qq, then it is preserved under evolution to larger 
values of |117l [T5] , but care must be taken since there are fits to the unpolarized data with 
parton densities which are negative at = IGeV/c^ [118j ! This is possible because in the 
unpolarized case one has the luxury of using only data typically with > 4GeV/c^, by which 
point the density has evolved into a positive one. 

• There is some guidance as to the behavior as x — 1. Perturbative QCD arguments |112j suggest 
that quark densities vanish like a power of (1 — x) which depends on the helicity; see the discussion 
in Section fTM In practice the exponents of {1 — x) are left as free parameters to be determined 
in the fit to the data, and tend to come out close to the expected values. 

• What the behavior should be as x — is not very clear. Based on the correspondence between 
X — >■ in the structure functions, and energy i/ — > oo in virtual-photon Compton scattering, one 
would expect the power of x to be controlled by the leading Regge intercept i.e 

(7i(x)~x-" (68) 

where a is the intercept of the ai(1260) meson Regge trajectory in the isovector channel, and the 
/i(1285) meson trajectory in the isoscalar channel [119j . These intercepts are not well known (see 
discussion in [13] ). Roughly 

-0.4 < aa, ^ Of, < -0.18 . 

However, we know that for unpolarized DIS the effective power a increases with for the range of 
X explored thus far, contrary to the Regge picture, and this is explained as due to QCD evolution 
(soft Pomeron vs. hard Pomeron). Bass [120] suggests that something similar might be expected 
for the polarized case. 

Study of the small x behavior of the DGLAP evolution equations |121j yields parton densities that 
grow more rapidly than any power of Iu^xq/x) (where Xq is the upper limit where the treatment is 
valid), and the growth rate increases with Q^. On the other hand summation of double logarithmic 
terms at small x, which are not included in DGLAP, leads to the expectation Cs x~"^ and C„s x"""' 
for the singlet and non-singlet contributions to gi respectively (The coefficients C are sensitive to 
the structure of the parton densities). This was originally studied in |122l 1123] using a constant, 
non-running value of and later by Ermolaev, Greco and Troyan |124j using a running coupling. 
Both approaches give the same structure with > a„s, but with somewhat different numerical 
values of the powers of x. In |124] the intercepts are 

as ^2ans^ 0.86, 

implying that the singlet flavor combination {gi + gi) should diverge more rapidly than the non- 
singlet combination [gi — g^) as x goes to zero. Thus one would expect 

+ > as x->0. (69) 

The only data on both flavor combinations at extremely small x values and > IGeV^ are 
from the SMC experiment |125] shown in Fig. [231 where precisely the opposite trend is visible. 
COMPASS data on the deuteron further conflrm the smallness of the flavor singlet (see Fig. [TT]) . 
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Further information will become available when COMPASS publish their new data on g^. How- 
ever, as experiments probe smaller and smaller x, there ought to be a dramatic change in the trend 
of the data, or else the theoretical arguments are incomplete. Indeed Blumlein and Vogt |126] 
suggest that corrections for sub-leading terms could significantly change these results. (A more 
detailed discussion can be found in Section 11.8 of [93]; for possible connections to the infrared 
renormalon model see the work of Kataev [127j ). 
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It is also argued |110j that 
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(70) 



Again, in practice, the power of x is taken as a parameter to be determined by the fit to the data. 

As explained in Section 2.2 there is, alas, no absolute physical meaning to parton densities in 
QCD. One has to specify in what factorization scheme they are determined. Although the differ- 
ences between densities in different schemes are proportional to ^^(Q^), these differences can be 
surprisingly large for the smaller densities e.g. As(x)|;^ ^ 2 ^s{x)\jet at = IGeV/c^ 



Although there has been a tremendous increase in the amount and quality of polarized data in 
the past few years, it is still usual to impose, rather than test, the various sum rules (Bjorken, 
SU (3) flavour) which should be satisfied by the quark densities. For a discussion of sum rule tests, 
see Section [321 

It is important to appreciate that the QCD formula for gi{x, Q"^), Eq. fl58l) . is a leading twist{LT) 
expression, i.e it neglects terms of order M'^/Q'^. For experiments at relatively low values of 
the expression Eq. (1S5]) must be supplemented by higher twist (HT) terms. There are two types 
of HT terms: target mass corrections, which are purely kinematical, and dynamical corrections. 
The target mass corrections can be calculated exactly and were first given in a closed form for gi 
by Piccione and Ridolfi |128j . For a new approach, see also |129j . The details and evolution 
of the dynamical HT terms are not known, so they are simply parametrized in the form h{x)/Q'^, 
with h{x) determined, after applying the target mass corrections, from experiment at a few values 
of X. For access to the literature see |130j . 

It is well known that parton densities extracted in LO and NLO sometimes differ significantly, 
and one may therefore worry about the effect of NNLO etc. corrections. Some indication of the 
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reliability of the NLO approximation can be obtained via the LSS Transformation Test [91]. In 
each order of perturbative QCD there exist transformation formulae relating parton densities in 
different factorization schemes. Let us indicate this symbolically for two schemes, A and B: 

Aq{x)\B = TB^A^qix)\A. (71) 

Suppose now that Tb^a is known to NLO accuracy, and the parton densities are extracted from 
the data, independently, in NLO, using schemes A and B, with results Aq{x)\'^^B 5 respectively. If 
the densities can be extracted reliably in NLO i.e NNLO effects are unimportant, then one should 
find 

Aq{x)\^^'- = Tb^a Ag(x)|t*^ (72) 
Any failure of this equality is a measure of the importance of NNLO effects. Thus the ratio 

Ag(x)|g'*'' - Tb^a Ag(x)|^°'° 
A(7(x) I ^"t" + Ts_ A Ag(x) I 

gives some indication of the reliability of the parton densities. 
A detailed discussion of many of the above issues can be found in |131j . 



2) Extraction of parton densities from DIS 



Following upon the remarkable EMC experiment in 1988, there have been many new data on po- 
larized DIS and many QCD analyses. We shall only be concerned with the recent analyses which 
have included all or most of the present world data, namely. Leader, Sidorov and Stamenov (LSS'05 
and LSS'06) [7l[l30], Alexakhin et al (COMPASS'06) [19], Hirai, Kumano and Saito (AAC'03) [T32] . 
Bliimlein and Bottcher (BE) |133j and Gliick, Reya, Stratmann and Vogelsang (GRSV2000) [108j . 
We shall also comment on the ground-breaking paper of de Florian, Sassot, Stratmann and Vogelsang 
(DSSV) [1U4] which analyzed simultaneously data on DIS, SIDIS and pion production in polarized pp 
collisions at RHIC, and on the work of Hirai and Kumano (AAC'08) [134] who analyzed data on DIS 
and on pion production at RHIC. Details of the parametrizations etc. can be found in the papers cited 
and also on the Durham HEP Databases website [135j . 

One of the problems in any analysis is to decide whether to use the lowest twist (LT) expression 
Eq. (135]) or to include higher twist (HT) terms. This is not just a question of the kinematic region in a 
given experiment, but also of what quantity one fits. Thus long ago we discovered empirically that HT 
terms cancelled out in the ratio 1^. Put 

9i - 9i +9i ^1 -^1 +^1 

Then 

r,LT „HT rpHT 

rgi_iEXP _ 9i_r. 9l fj_l 

L^^J - pLTV'+ gLT pLT\ 

(73) 



9^^ 



pLT 

„HT pHT 

provided there is a cancellation between ^^xr and -^ttt- 

The left panel of Fig. [21] which compares HT terms for gi [130j and for Fi [136] demonstrates the 
validity of this for x > 0.15, but clearly indicates that ignoring HT terms in the ratio |r below x = 0.15 
is incorrect. The HT terms for gi determined in [130j are shown in the right panel of Fig. [2l] using 5 
bins in x. (The results change negligibly if 7 bins are utilized.) 
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Figure 24: Comparison of HT terms in gi and Fi (left panel) and higher twist terms for 
protons and neutrons extracted from data (right panel). 



A more subtle problem relates to the kinematic relations between various quantities. The exact 
relations between A\\, Ai and gi include HT terms involving 7^, both in Eq. fl52]) and in the expression 
for D, Eq. fl2Bl) . If it is safe to ignore HT terms in gi, Fi is it safe to ignore them in these kinematic 
relations as well? DSSV assume it is, whereas LSS insist that if 7^ is significant such terms must be 
kept in the kinematic relations. The validity of the latter viewpoint is borne out by the empirical fact 
(accidental?) that HT terms can be ignored in |r- in the region x > 0.15, where, in some cases, 7^ is 
non- negligible. Thus the DSSV analysis, which treats |r- in LT, but ignores 7^, may be unreliable in 
those regions where 7^ is important. This issue is particularly important for data from HERMES and 
JLab. 

• The light quark densities. There is broad agreement between the various analyses for the Au{x) + 
Au{x) and Ad{x) + Ad{x) parton densities. Generally the constraint Eq. ( l66l) has not been 
enforced in the parametrizations used for the d quark . Early data demanded negative values of 
Ad{x) + Ad{x) and continued to do so even when the measured region was extended to x = 0.6 
at Jefferson Laboratory [611 [72] . These data are discussed in detail in Section 12.41 Note that 
a fit [137j using the Brodsky, Burkardt and Schmidt type parametrization |110j . which respects 
Eq. (1^ . led to a Ad{x) + Ad{x) which became positive just beyond x = 0.6. With the 12GeV 
upgrade at Jefferson it should be possible to explore out to x = 0.8 and to settle the matter. Fig. [231 
shows a comparison between the LSS'06, DSSV and AAC'08 Au + Au and Ad + Ad densities. 
It is seen that the agreement is impressive, although the d-quark polarization Ad{x)/d{x) is still 
rather uncertain at high x (see also Section [274]) . 

In Fig. [26l we show the impact of the CLAS'06 and COMPASS'06 data on the uncertainties in 
x{Au + Am) and x{Ad + Ad) which had been found in the LSS'05 analysis. The uncertainties are 
estimated by the usual Hessian method, utilizing Ax^ = 1. It should be borne in mind that these 
uncertainties, and all those to follow, must be treated with caution. They are calculated relative 
to a chosen analytical form for the parametrization of the parton density and do not at all reflect 
possible uncertainties in the form of parametrization. 
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Figure 26: Impact of CLAS'06 and COMPASS'06 data on the uncertainties in x{Au + Am) 
and x{Ad + Ad) of the LSS'05 analysis. 



• The polarized strange quark density. This is a highly controversial issue at present. All analyses 
of purely DIS data have found negative values for As{x) + As{x). Fig. [771 shows the impact of the 
CLAS'06 and COMPASS'06 data on the uncertainty in the LSS'05 determination of x{As{x) + 
As{x)). 

An important quantity is the first moment 

AS= [ dx[As{x) + As{x)]. (74) 
Jo 

LSS'06 [no] give for its value 

A5^= -0.126 ±0.010 at Q^ = lGeV'^ (75) 

and in |138j it was shown that a positive value for the first moment would imply a huge breaking of 
SU{3)f invariance, far greater than the ±10% breaking estimated by Ratcliffe |139] . Nonetheless 
analyses of SIDIS data for kaon production have suggested positive values of As{x) + As{x) 
for X > 0.03. A hint of this first emerged in the HERMES analysis of 2004/5 [140j and was 
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Figure 27: Impact of CLAS'06 and COMPASS'06 data on uncertainties in the LSS'05 de- 
termination of x(As(x) + As(x)). For curve labels see Fig. [261 

confirmed by the more precise data reported in [102] where the first moment for the measured 
range 0.02 < a; < 0.6 was given as 

AS = 0.037 ± 0.019(stat) ± 0.027(sys.) (76) 

One could argue that the HERMES results should not be taken too seriously, since the analyses 
were carried out in LO and relied on derived quantities called purities, whose accuracy may have 
been overestimated. However, the recent combined analysis of DIS, SIDIS and pp 7r° or jet + X 
by the DSSV group [104j . which is an MS NLO, LT analysis, also finds positive values for 
As(x) + As(x) for x > 0.03, yet ends up with a negative first moment AS" = —0.114 at = 
lOGeV"^. This is achieved by As(x) + As{x) becoming negative below roughly x = 0.02. But 
there are essentially no data in the latter region, which suggests this must be caused by the need 
to satisfy the SU{3)f symmetry condition Eq. (H^ . DSSV state that they do not impose SU{3)f 
symmetry, and they multiply the RHS of Eq. (H6|l by (1 + esu{3)) where esu{3) is a free parameter. 
However, a bias against large values of esu(3) is built into the minimization procedure which 
explains the turning negative of As{x) + As{x), although this is not stated in the quoted paper. 

The recent AAC'08 analysis [134j . based on the world data on DIS plus the vr" production data 
from RHIC (incorporated at LO), also finds a negative strange quark density. Fig. [251 compares 
the LSS'06 version of As + As with that of DSSV and AAC'08. Note that the shape appears 
somewhat different in the various DIS analyses. The origin of this has been traced to the different 
choices of unpolarized densities used in the imposition of positivity |141j . as can be seen in the 
figure. 

It is thus clear that the origin of the positive values for As{x) + As{x) lies entirely in the SIDIS 
data. It may be significant that both the HERMES and DSSV analyses utilize the recent de 
Florian, Sassot and Stratmann [142j set of fragmentation functions, which are surprisingly different 
from earlier sets of fragmentation functions, especially for kaons. COMPASS are studying the 
dependence of As{x) + As{x) on the choice of fragmentation functions. A preliminary report can 
be found in the talk of R. Windmolders at the SPIN2008 Conference. 

Finally, there is an intriguing question as to whether the standard coUinear analysis should be 
applied to the SIDIS data. The point is the following. What we mean by any q{x, Q"^) or Ag(x, Q"^) 
actually involves an integral over the intrinsic transverse momentum kt of the quark up to a 

private communication from Werner Vogelsang 
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Figure 28: LSS'06, DSSV and AAC'08 (MS) versions of x(As + As). 



maximum k"^ ^ Q"^. If, as is the case for some of the SIDIS data, the magnitude of the transverse 
momentum of the detected hadron is less than Q, then it is perhaps incorrect to analyze the data 
using the conventional coUinear densities. Indeed, Bass [143j has shown that for that part of the 
strange quark density generated via photon-gluon fusion the sign of AS" is sensitive to the upper 
limit in the integration, but the overall effect is much too small to be relevant. 

• The flavor singlet first moment AS. All the modern global analyses obtain compatible values for 
AE. In the WS scheme, where aoiQ'^) = AS(Q2) they find at = A GeV^: 



LSS'06 


COMPASS'06 


AAC'08 


DSSV 


0.24 ±0.04 


0.29 ±0.01 


0.25 ±0.05 


0.24 



For reasons which are not understood these values are somewhat lower than the values obtained 
directly from Tf as discussed in Section 13.21 

• The polarized gluon density. In Sec. 12.21 we explained how the "spin crisis in the parton model" 
could be alleviated by generating the experimentally small value of ao through a cancellation 
between relatively large values of AS and AG (see Eq. I62l) . Present day estimates are oq ~ 0.2. 
Thus demanding AS ^ 0.6 requires, for the first moment, 

AG ^1.7 at Q'^ = lGeV'^ (77) 

The question is whether this is compatible with what we know about the polarized gluon density. 

There are three ways to access AG(a;): via polarized DIS, and via the measurement of the asym- 
metry All both in SIDIS production of charmed quarks or high px jets, and in semi-inclusive 
polarized pp reactions at RHIC. Here we will only discuss DIS. The latter approaches have been 
treated in Section 12.31 

In determining the parameters of the polarized gluon density from fits to the data on gi{x, Q^), it 
should be remembered, as already mentioned, that the main role of the gluon is in the evolution 
with Q'^, but that the range of is very limited, so the determination of AG(a;) is imprecise. 
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Figure 30: Left: a compilation of earlier results for xAG{x)] right: present situation. 

2 

For a long time all analyses seemed to indicate that AG{x) was a positive function of x. LSS used 
a very simple parametrization 

xAG{x) = ^AgX^'lxGix)] (78) 

where Ag is fixed so that the first moment of AG is given by r]g. In the minimization procedure 
there was nothing to stop r]g from being negative. Yet the best values always corresponded 
to positive AG{x). With the inclusion of recent data, LSS'06 |13U] find equally good fits with 
positive, negative and sign-changing AG{x), provided higher twist terms are included. The latter 
are particularly demanded by the CLAS'06 data [72]. Fig. [29] compares the LSS'06 positive, 
negative and sign-changing versions of AG{x). 

Also the COMPASS'06 analysis ^9], even though HT terms are not included, finds acceptable 
negative AG{x) fits, but has some peculiarities (linked to the fact that the flavor combinations 
Ag3 8 and AS are parametrized, rather than the individual quark densities), which suggest it may 
not be very physical. The DSSV analysis finds a sign-changing AG{x), while the AAC'08 analysis 
finds both positive and sign-changing densities. In the left hand panel of Fig. [30] we show various 
results for AG{x) as of a couple of years ago. The right panel shows the present situation. In 
Fig. [31] we show the impact of the CLAS'06 and COMPASS'06 data on the uncertainty in the 
LSS'05 determination of xAG{x). 
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Figure 31: Impact of CLAS'06 and COMPASS'06 data on uncertainties in the LSS'05 de- 
termination of xAG{x). For curve labels see Fig. [261 



Regarding the present data, in all fits, and irrespective of the form of the gluon density, the 
magnitude is always found to be very small. Typically one has | AG| ^ 0.29 ± 0.32, much smaller 
than the desired 1.7 ! As discussed in Section [231 ciU the present data on SIDIS and pp reactions 
are perfectly compatible with this very small magnitude for AG and cannot distinguish between 
the various sign possibilities. 

2.6 The Spin Structure Function g2 

The spin structure function g2, unlike gi and Fi, does not have an intuitive interpretation in the simple 
quark-parton model. To understand g2 properly, it is best to start with the operator product expansion 
(OPE) method (see also Sections 12. 1112. 5113. 1113. 41) . In the OPE, neglecting quark masses, g2 can be 
cleanly separated into a twist-2 and a higher twist term: 



(this holds also with the inclusion of TM corrections |144j ). and the higher-twist term g2''^'{x, Q"^) arises 
from quark-gluon correlations. Therefore, g2 provides a clean way to study higher-twist effects. 

Experimentally, g2 can be extracted from combined measurements of both the longitudinal and 
the transverse target spin cross section differences (Eqs. [T4][T5] in Section II. ip . or the corresponding 
asymmetries A\\, A^. These measurements demand high precision and large luminosity, since the factors 
multiplying g2 tend to be relatively small in DIS kinematics. Several experiments pT| I145[ 1146] 
therefore took only enough statistics on the transverse target configuration to reduce the systematic 
error on gi due to the unknown contribution from g2 below the statistical error, yielding only crude 
information on g2 itself. 

Two experiments at SLAG, E143 and E155 [5^ W\\ took dedicated data to measure (?2 and the 
virtual photon asymmetry A2 (see Eq. [23]) with more precision. The most extensive data set in the 
DIS region comes from the second run of E155, both for the proton and the deuteron, with g2 for the 
neutron extracted from the difference. The general trend of these data follows the Wandzura-Wilczek 



g.M')=g7''M')+9rM')- 



(79) 



The leading-twist term can be determined from gi as [H] 




(80) 
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Figure 32: Results for g2 of the proton and the deuteron (top-left panel) and the neutron 
(top-right panel) as a function of x from SLAG ^1] and JLab Hall A [621 163] , in comparison 
with (solid curves). The left panel also shows bag model calculations of Stratmann |147] 
(dash- dot- dot) and Song [148j (dot) and the chiral soliton models of Weigel |149[ 1150] (dash 
dot) and Wakamatsu |151j (dash). The bottom plot shows results for g"^ from JLab Hall 
A [S3] as a function of at an x value around 0.2, in comparison with g^^ and models. 



form (Eq. [HOl) . but some deviations are seen for the proton, especially at small x. The data on A2 show 
that this asymmetry is indeed rather small, well below the limit given in Eq. [301 

Figure [321 shows the SLAG results on the proton, deuteron (top-left panel) and neutron (top-right 
panel) as a function of x. Also shown on the top-right panel of Fig. [321 are recent JLab neutron 
results with polarized ^He from Hall A E99-117 [62] and E97-103 [63]. E97-103 is a dedicated precision 
measurement of g2 in the DIS region at low to study its dependence. It covered five different 
values from 0.58 to 1.36 GeV^ at x ~ 0.2. Results on the dependence of (?2 ^^^^ given in the bottom 
panel of Fig. [321 The light-shaded area in the plot gives the leading-twist contribution, obtained by 
fitting world data |133j and evolving to the values of this experiment. The systematic errors are 
shown as the dark-shaded area near the horizontal axis. The precision reached is more than an order 
of magnitude improvement over that of the SLAG data [H] . The difference of g2 from the leading twist 
part (g^^) is due to higher-twist effects and is sensitive to quark-gluon correlations. The measured g2 
values are consistently higher than g^^. For the first time, there is a clear indication that higher-twist 
effects become significantly positive at below 1 GeV^, while the bag model |147] and Ghiral Soliton 
model [1501 1151j predictions of higher-twist effects are negative or close to zero. The g^ data obtained 
from the same experiment agree with the leading-twist calculations within the uncertainties. 
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2.7 Spin Structure in the Resonance Region 



While most measurements of spin structure functions initially focused on the deep-inelastic region, 
more recently experiments have also collected a truly massive data set on gi and g2 in the region W < 2 
GeV, the nucleon resonance region, over a wide range in Q^. Because the (resonant) final state in this 
kinematic region can have a strong influence on scattering matrix elements, these data are less suitable 
for a direct interpretation within the framework of perturbative QCD (Section 12.21) . On the other hand, 
they can reveal new and important information on the internal structure of nucleon resonances and 
their excitation by the electromagnetic probe. Combining proton, deuteron and neutron targets, one 
can separate the various spin-isospin channels that lead to the excitation of overlapping resonances as 
well as non-resonant (multi-) meson production. 
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Figure 33: Data on the spin structure functions gi (filled circles) and g2 (open squares) vs. 
X in the resonance region, for ^He (approximately a neutron target), from Jefferson Lab's 
Hall A. Statistical errors are shown as error bars, while the systematic errors are indicated 
by the shaded bands. The average for each data set is indicated in the figure. 



For a given resonance, the (electro- or) photo-excitation strength can, in general, be given in terms 
of three helicity amplitudes, namely A-^/2{Q'^) (transverse photons leading to final state helicity 3/2), 
^1/2 (Q^) (transverse photons leading to final state helicity 1/2) and 5*1/2 (Q^) (longitudinal photons). 
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These amplitudes are directly related to the photon asymmetries, vz. 

Ai - 



\Al/2\ 


2 _ 


1^3/2 
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1^1/2! 
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^^3/2 


2 



. ^ '^1*72-41/2 

^ g*|^i/2P + l^3/2r 

Here, q* is the (virtual) photon three-momentum in the rest frame of the resonance. As a consequence, 
electromagnetic excitation of spin-1/2 resonances will have an asymmetry Ai = 1 since the amplitude 
A3/2 cannot contribute. Vice versa, for resonances like the Delta (with spin 3/2) that are predominantly 
excited via Ml transitions, one has ~ -\/3^i/2 and therefore Ai ^ —0.5. By studying the Q^- 
dependence of these asymmetries for a given mass range, one can gain information on the relative 
strength of overlapping resonances (and non-resonant background) as well as, for spin-3/2 resonances, 
the transition from A3/2 dominance at low to the asymptotic A1/2 dominance expected from helicity 
conservation in coUinear pQCD. 

Measurements of spin structure functions in the resonance region are also required for the evaluations 
of moments (Section [3]), since in most cases the integrals must be evaluated up to x = 1, which includes 
the resonance and elastic region, particularly at more moderate Q^. As an example, for = 3 GeV^, 
half of the x-range, 0.5 < x < 1, lies in the resonance region. These measurements are also very 
important to study the onset of duality for spin structure functions, see Section HI Finally, higher order 
QED radiative corrections to all measurements (even in the DIS region) have contributions from the 
resonance region, making good knowledge of spin structure functions at low W imperative for accurate 
interpretation. 

The first experiment to measure double spin observables at SLAG was also the first to take data in 
the resonance region [2U]. The next dedicated measurement was by the E143 collaboration [32], who 
used a lower beam energy of 9.7 GeV to cover the resonance region at = 0.5 GeV^ and 1.2 GeV^. 
These data showed that the asymmetry Ai is indeed negative in the A region at moderate Q^, while 
it rises rapidly to rather large values in the second resonance region (dominated by the Sn and D13 
resonances) . 

Since then, experiments in all three halls at Jefferson Lab have collected a vast data set on all 
inclusive spin structure functions throughout the resonance region, for proton, deuteron and ^He targets, 
spanning several orders of magnitude in (from 0.015 to 5 GeV^). Figures 1551 through [5^ show some 
samples from this data set. Double spin asymmetries in the resonance region were also measured at 
MIT Bates with the BLAST detector |152j on the proton and the deuteron at 0.85 GeV beam energy. 
The kinematic region covered ranges from = 0.08 — 0.38 GeV^ and from the elastic peak toW = 1.33 
GeV. More data in the region of the A(1232) resonance were also obtained at NIKHEF for = 0.11 
GeV2 [T^. 

The data in Fig. [33] are from experiment E94-010 in Jefferson Lab's Hall A [57]. Both spin struc- 
ture functions gi and g2 were extracted directly from cross section differences for polarized electrons 
scattering off a ^He target polarized both along and transverse to the beam (see Eqs. [MfTSl) . These 
data approximate the spin structure functions for a free neutron, although nuclear corrections were not 
applied. One can clearly see the strongly negative response of gi in the A (1232) resonance region (high 
x). The structure function g2 is close to the negative of gi, which implies that A2, which is proportional 
to gi + g2, is close to zero. 

Fi gure [Ml shows a sample from the vast data set collected with GLAS during the EGl experiment [72] 
on the spin structure function gi for the proton and the deuteron. gi was extracted from measurements 
of the asymmetry A\\ (see Eq. [HT]) over a range in x covering both the resonance and deep inelastic 
regions, for over 25 bins in Q^. As increases over 2 orders in magnitude, the asymmetry (and 
therefore gi) changes from a predominantly negative value (/13/2-dominance) to the strongly positive 
DIS limit for nearly all x except in the region of the A(1232) resonance. 
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Figure 34: Data on the spin structure functions gi of the proton vs. x in the resonance 
region and beyond, from Jefferson Lab's Hall B. Statistical errors are shown as error bars, 
while the systematic errors are indicated by the shaded bands. The average for each data 
set is indicated in the figure. The curves are from a phenomenological parametrization of 
the data. A similarly complete data set exists also for the deuteron. 



Finally, we show in Fig. [35] the data from the "RSS" experiment in Hall C [HQ]. This experiment 
measured gi and g2 in the resonance region for one bin in Q^, centered around 1.3 GeV^, for both 
proton and deuteron, using transversely and longitudinally polarized targets. The gi data "oscillate" 
around the deep inelastic curves (from various NLO fits), while g2 is in rough agreement with the 
Wandzura-Wilczek form (not shown). 
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Figure 35: Data on the spin structure functions gi and g2 of the proton versus x in the 
resonance region, for an average of 1.3 GeV^, from Jefferson Lab's Hall C (dark filled 
symbols) and other experiments (as indicated). Statistical errors are shown as error bars, 
while the systematic errors are indicated by the shaded bands. 



3 Sum Rules and Spin Polarizabilities 

3.1 Moments of Spin Structure Functions 

Moments of structure functions are a powerful tool to study fundamental properties of nucleon structure, 
like the total momentum fraction carried by quarks or the total contribution of quark helicities to the 
spin of the nucleon. While a complete description of structure functions based on fundamental QCD 
principles may be unattainable for now, moments of structure functions can be directly compared to 
rigorous theoretical results, like sum rules, lattice QCD calculations and chiral perturbation theory. 
The original "spin crisis" was directly due to a discrepancy between the data and an approximate "sum 
rule" by Ellis and Jaffe [3]. On the other hand, a detailed measurement of the Q^-evolution of the 
Bjorken sum rule [5] provides a significant test of pQCD in the spin sector. 

As mentioned in Section 12. via the Operator Product Expansion (OPE) the moments of gi^2 can 
be related to hadronic matrix elements of current operators. The moments are given as a sum, ordered 
according to the twist r = (dimension — spin) of the current operators, beginning with the lowest twist 
r = 2. Each additional unit of r produces a factor of order ^"^^ , and is thus less important in the high 
q2 j-QgiYYie. Strictly speaking, the higher twist terms are mixed up with correction terms of order 
of purely kinematic origin (target mass corrections), which can be calculated exactly, so as to expose 
the HT terms of dynamic origin. Some of the twist 2 terms are directly measurable in other processes. 
Twist 3 and higher terms can sometimes be determined from combinations of measured quantitiesJl 

After a long series of manipulations one finds expressions for the n*'^ moment dx x"^'^ gi{x,Q'^) 

for n = 1,3,5 of gi and the n^^ moment dx x'^~^ g2{x,Q'^) for n = 3,5,7 of g2, in terms of 

hadronic matrix elements of certain local operators. (For details see Section 22.2 of |156j ). 

The most important case is the first moment, n = 1, of gi, where the operators involved are the 



^The reader is warned that the notation in the Uterature on this subject is in a state of anarchy, with identical symbols 
often being used for significantly different quantities in different papers. Thus Gi^2 sometimes means times our G1.2, 
but perhaps more bizarre is the use of the symbol for P • g in [1541 1155j 
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octet of well known axial-vector currents which control neutron and hyperon /5-decay: 

Jl, = i^w(^Y)'^ (^ = l,2,...,8), (82) 
where the Xj are the usual Gell-Mann matrices and is a column vector in flavor space 

= MAd • (83) 




The flavor singlet current, which does not play a role in /9-decay, is 

= V'7m75V' • (84) 

The relation of the hadronic matrix elements to the first moments = qa, dg, clq of the flavor 
combinations of quark densities introduced in Eqs. HH is 

{P,S\Jl^\P,S) = Ma3S, 
{P, S\J^^\P, S) = asSf, 
{P, S\Jl\P, S) = 2MaomS,. (85) 

The dependence of Cq arises because the singlet current has to be renormalized and it is customary 
and convenient to choose as the renormalization scale. 

At leading twist the expression for Ff'*^, valid for ^ M^, can be written 

rr = ^ [( ± as + ^ as)E^s{Q') + ^ ao{Q')Es{Q')] , (86) 
where the coefficient functions are given, for 3 active flavors, by 

Ens{Q') = !-(-)- 3.58 {^f - 20.22(^)^.. (87) 



and in the MS scheme 

= !-(-)- 1.096 ( 

71 71 

Note that the singlet coefficient here multiplies ao(Q^) and not the invariant quantity referred to in 
various papers as cio, g^^linv and Al^linvariant, which are the value of ao(<5^) as oo. 
There exist important sum rules for the moments of gi^2 which we now discuss. 

• The Bjorken sum rule. Using only isospin invariance Bjorken [5] showed that as oo, 

/' dx [^?(x, Q') - g^{x, Q')] = ^ [1 - (^) - 3.58(^)^ - 20.22(^)^ ■ ■ .]., 
Jq d tt tt tt 



^9) 



where the square bracket on the RHS contains the perturbative corrections Eq. [871 for 3 active 
flavors, to Bjorken's original result qa/Q- The sum rule is quite rigorous, involving only the 
assumption of isospin invariance, and inspired Feynman to say that its failure would signal the 
demise of QCD. 
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The Efremov, Leader, Teryaev sum rule . The ELT sum rule [TT] is also a rigorous result, following 
from the OPE, and states that 



[ dxx[glf{x) + 2glf{x)]=0, 
Jo 



(90) 



where V implies valence contribution and the result holds for each flavor /. It was originally 
thought reasonable to assume that the sea-quark densities are the same in protons and neutrons, 
in which case Eq. [90] can be written as a kind of analogue to the Bjorken sum rule 

' dxx [g',{x) - g^{x)] = \ j'^ dxxig'lix) - gl{x)]. (91) 

However, the above assumption is tantamount to assuming Am = Ac/ in a proton, which, at least 
for the unpolarized antiquark densities, is not a good approximation. For further interesting sum 
rules of this type see [144^ 1157] . 



The Burkhardt-Cottingham sura rule. The BC sum rule [10] 

dxg2{x) = Q (92) 





does not follow from the OPE, although it looks as if it does, because the n^^ moment of g2 is 
proportional to (ra — 1); but, as mentioned, the OPE result for g2 only holds for > 3. 

The BC sum rule is based on an assumption about the asymptotic behavior of a particular virtual 
Compton amplitude (see discussion below). The argument is not watertight and it may be that 

g2{x) as X 0, (93) 

thus making the integral in Eq. [92] diverge. For a detailed discussion of this issue, see Section 5.2 
of [T2|. Note that the Wandzura-Wilczek approximation to g2, Eq. [38] automatically satisfies the 
BC sum rule 

1 

dxgf'^{x)=Q. (94) 







Higher twist corrections. 

At lower values of higher twist contributions become important. Of particular interest is 
ri((5^), which can be expanded, for Q"^ > Ag^^, in inverse powers of (twist expansion) 

Tm = ^M')+f^+^+ (95) 



with /i2, the leading twist term, given by Eq. 

The evolution in of the ^niQ"^) can, in principle, be calculated in perturbative QCD, but 
results are known only for /i2 , and it is usual, therefore, to ignore the dependence of fi^^Q. 
The functions fin are related to matrix elements of operators of twist t < n; the occurrence of 
operators with twist lower than n is a kinematic effect, a consequence of target mass corrections. 

We consider now the structure of the higher twist (i.e., 1/Q^) corrections to Ti^Q'^) as defined in 
Eq. [95] As mentioned [JiiiQ'^) contains both target mass (TM) corrections and genuine dynamical 
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higher twist contributions. The TM corrections for (71,2(2;, Q^) were given in [128] . For the first 
moment they yield 



(96) 



g2, as discussed in Section 12. 6[ contains a twist-2 (and kinematic twist-3) part, , given in 
Eq. [8OI and a dynamic twist-3 part given by g2 — g^^ ■ Writing 

92 = 97"^ + {92 - 9l^) (97) 

one finds, after a httle algebra, that 



)] ■ (9J 



^T\Q') = ^{ rfxx^ ^i(x, Q2) + 6 dxx' [^2(x, g^) - ^7r^(x, (f 

The twist-2 combination of moments is usually referred to as 02 and the twist-3 as d2. More 
precisely 

a2{Q^) = 2 [ dxx^gi{x,Q^) d2{Q^) = 3 [ dxx^g2{x,Q^) - g^'^ix^Q^ 

Jo Jo 

There is also a genuine twist-4 contribution to fx^, written as ^^/2, so that finally 



(99) 



fi,{Q') = —[a2{Q') + 4d2{Q') + Af2{Q% (100) 

The twist-4 part /2 cannot be written in terms of moments of the standard scaling functions. In 
terms of operator matrix elements it is defined by 

f2{Q')M'S^ = ^ J2 e}{P,S\g^PfG^''^.^|Jf\P,S) (101) 

flavors 

where G^'^ = [1 / 2)e^'^"^ G with the sign convention e^^^^ = -|-1. ( Note that in some papers the 
normalization of the covariant spin vector is taken to be 5*^ = — M^, in which case will appear 
on the LHS of Eq.fMi) 

The first correct expressions for d2 and /2 were given by Ji and Unrau |155j and by Ehrnsperger, 
Mankiewicz and Schafer [158] . ( Note that in Ji's paper the twist-4 term in Eq. llOOl appears with 
a minus sign. This is because his coupling g in the operator in Eq. llOll has the opposite sign to the 
usual convention i.e. he uses for the covariant derivative acting on quark fields = + igA^.) 

Since /i4, 02 and ^2 can be measured, Eq. llOOl yields a measurement of the twist-4 matrix element 
f2- Both d2 and /2 can then be compared to non-perturbative QCD calculations. 

Finally, it is of some interest to consider the color electric and color magnetic polarizabilities of 
the nucleon |159j . They are defined by 

Xe2M^S ={N,S\jaXE, \N, S ) Xb^M^S ={N,S\ j^B^ |iV, S ) (102) 

where S is the rest frame spin vector of the nucleon and is the quark current. Ea and Ba are 
the color electric and magnetic fields respectively. The color polarizabilities can be expressed in 
terms of ^2 and /2 as 

XE=l{2d2 + f2) XM = \{^d2- f2). (103) 
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• Extrapolation to low Q^: the Gerasimov, Drell, Hearn sum rule. The twist expansion does not 
converge for very small Q^. Thus to study polarized structure functions in the low region it is 
necessary to make contact with soft, non-perturbative physics i.e. to study the relation to Comp- 
ton scattering with real, or almost real photons. In this region care must be exercised in defining 
moments of structure functions. As already mentioned in Section 12.11 relations between moments 
of structure functions and matrix elements of operators are only valid if the moments include the 
elastic contributions located at a; = 1. In the deep inelastic region the elastic contributions are 
negligible and are not included in experimental estimates of the moments, but at low Q^, in the 
resonance region, the distinction is important. Thus moments in the latter region, without an 
elastic contribution will be labeled T. 

The GDH sum rule |160l 1161] concerns Compton scattering with real photons, i.e., with Q"^ = 
0, and follows from the fact that the value of the forward spin-flip amplitude, f2{i^), a.t u = 
0, calculated to order e^, is given by low energy theorems, and from the assumption that the 
dispersion relation for it does not require subtractions. This leads to 



where k is the anomalous magnetic moment of the nucleon. The aA,p are the total cross-sections for 
the absorption of a circularly polarized photon by a proton polarized with spin antiparallel/parallel 
to the photon spin. 

Since, to order e^, the cross-sections are zero below pion production threshold uq (see below), and 
following the convention used by experimentalists, we write Eq. 11041 in the form 



In the language of DIS the above cross-section difference is referred to as arr and Eq. I105[ 
generalized to non-zero Q^, is usually written 




(104) 




(105) 




(106) 



where 




(107) 




2M2 




with 




(109) 



and where TTriQ"^) is the inelastic portion of the first moment 





and Xo((5^) is the threshold for pion production. 




Xq = 



Q2 + m^(2M + m^)' 
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Note that the generahzation of /(O) in Eq. 11051 to arbitrary is then given by 



im = ^-^iTT{Q\ (112) 

The relevance of Eqs. 11061 and 11081 to deep inelastic scattering was first pointed out by Anselmino, 
loffe and Leader [162] . who noted that the first moment rj,j.((5^), which is a slowly (logarithmi- 
cally) varying, positive function in the DIS region, would have to change drastically in order to 
satisfy Eq. 11061 at small Q^. However |162j failed to distinguish between the full moment r^y((5^) 
and its inelastic version r^j,(Q^), and Ji |154j pointed out that there is a significant difference 
between the extrapolation to = of the full and the inelastic moments . The reason is a 
non-uniformity of the limits — > and — ^ in the generalization of f2{y) to ^ 0, namely, 

hm lim /2(i/,g2)^lim lim h{y.Q^). (113) 

Q2_>Q ;y_>0 i/— >0 Q2_,o 

Thus whereas the RH limit does not inherit an elastic contribution from the Born terms in 
Compton Scattering the LH one does. And since the integrals involved in the moments above 
correspond to z/ = in f^i^^jQ"^), it is the LH limit which must be used in the extrapolation to 
= 0. Thus for the full moment Eq. 11061 is altered to 

^hm IMQ') = Q^lFfiQ") + - '^V4 (114) 

where Ff^ are the Dirac and Pauli elastic from factors normalized to 

<(0) = 1 i^il(O) = Ff^^M = (115) 

In a similar way one can generalize to arbitrary |154j the spin-dependent Compton amplitudes 
Si^2, normalized so that 

ImSi,2{iy, Q') = 2n Gi,2{i^, Q^). (116) 
For 5*1 one obtains the dispersion relation 



Then 



-S,{0,Q')=h{Q^) = MM ^MG,{i^,Q^) (118) 



4 Jmq') ^ 

= ^r,(Q^) (119) 

where Ti^Q"^) is the inelastic portion of the first moment Ti^Q"^) 

ri(g2)= / dxg^{x,Q^) (120) 

^0 

and 

lim /i(Q^) = -K^/A. (121) 
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For the full moment, 



/l/f2 

lim h{Q') = -^Ftm [Ftm + F^\Q^)] - n'/A. (122) 

In the DIS region the elastic contribution to the moments is totally negligible. However in the 
extrapolation down through the resonance region towards = it is vitally important to dis- 
tinguish between the two cases. 

We turn now to the Compton amplitude 5*2(1^) which underlies the BC sum rule. If one assumes 
that 5*2(1^,(5^) satisfies a superconvergence relation, i.e., that it vanishes as z/ ^ cxd fast enough 
so that both 5*2 and VS2 satisfy unsubtracted dispersion relations, then one can show that 

/•oo /»oo 

/ ImS2{y,Q'^)du = 2ti j G2{iy,Q^)di^ = (123) 
Jo Jo 

which leads to Eq. [92] for the full moment. From a knowledge of the elastic terms in 6*2(0, Q'^) one 
obtains an expression for the inelastic integral 

hm ^Q^l 92{x, Q')dx = -Ff{Q')[Ff{Q') + F|'(Q2)]. (m) 

For a discussion of the convergence properties of the amplitudes see Section 2.10 of [163] . 

Generalization of the GDH approach. 

The above type of analysis, based on the dispersion relations for f2{ty,Q'^), Si^2{'^Q'^), has been 
generalized to all the amplitudes in virtual Compton scattering by Drechsel, Pasquini and Van- 
derhaeghen |164j . For very small energies, u < lyoiQ"^), the amplitudes can be expanded in powers 
of u"^. The coefficients of the next-to-leading terms are called generalized forward spin polarizabil- 
ities of the nucleon. They can be expressed in terms of moments of the structure functions and 
can thus be measured. They are important because they reflect soft, non-perturbative aspects 
of nucleon structure, and can be calculated, at least approximately, using various forms of chiral 
perturbation theory and lattice methods. They thus provide benchmark tests for these theories. 
Here we shall simply give the expressions for experimentally relevant polarizabilities. For details 
of the amplitudes etc, see [164] . 



00 



7o(Q^) ^ 7tt(Q') = 2a I ^ [uMG,{u, Q') - Q'G2{u, Q')] (125) 



a;2[^i(x,g2)____^2(x,Q')]rfx (126) 



Q 



6 



W)^^lt{Q'') = 2a r^[MG,{u,Q') + uG2{u,Q')] (127) 

J un ^ 

/ x'[gi{x,Q^)+g2{x,Q^)]dx (128) 



16aM2 /-^o 



The comparison between theory and experiment is discussed in the next few sections. 



48 



Table 1: Results for the first moments of the spin structure functions gi from different 
experiments. Each experiment evolved its data to a fixed value of which is indicated. 
The results for the deuteron are not corrected for its D-state, but are "per nucleon" . 



Q2 [GeV2] 


r? 


J- 1 


rf 


-pP pn 
J- 1 J- 1 


Ref. 


10 


0.120 ± 0.016 




0.019 ± 0.015 


0.198 ± 0.023 


SMC [2H] 


3 




-0.033 ± 0.011 






E142 [31] 


3 


0.133 ± 0.010 


-0.032 ± 0.018 


0.047 ± 0.007 


0.164 ± 0.023 


E143 [M] 


5 




-0 056 -1- 009 




168 + 010 


El 54^ l38l 


5 


0.118 ± 0.008 


-0.058 ± 0.009 




0.176 ± 0.008 


E1551 [3S] 


2.5 


0.120 ± 0.009 


-0.028 ± 0.009 


0.043 ± 0.004 


0.148 ± 0.017 


HERMES^ H] 


3 






0.046 ± 0.006 




COMPASS [19] 


5 


Bjorken Sum Rule 






0.182 ± 0.002 





^ From a NLO analysis 

^ Over measured region x > 0.021 only 



3.2 First Moment of gi at High and Sum Rules 

After the EMC announced their surprising result [2], the next round of experiments at SLAC, CERN 
and DESY was focused on the goal to measure the first moments of gi for the proton, deuteron and 
neutron (^He) in the DIS region with increasing precision. The main interest was to gain information 
on the contribution of quark helicities to the nucleon spin and to test the Bjorken sum rule [5]. 

Experimentally, these integrals can only be measured over some finite range in x, since the limit 
a; — > for finite requires infinitely high energy transfers u. In practice, all existing experiments 
accessing gi are fixed target experiments, with a lower limit of x of the order x = 0.004 — 0.01 (if 

> 1 GeV^). Experimental results for the moments Fi therefore include some extrapolation to x = 0, 
typically based on Regge phenomenology [1201 1165] I166j . NLO fits or some other assumption about the 
functional form of the underlying quark densities. This extrapolation to the lower integral limit is the 
largest source of theoretical uncertainty for the results obtained so far. Most experiments also use an 
extrapolation at high x, since small statistics and limited momentum resolution prevent an accurate 
measurement in this kinematic region. However, at sufficiently large Q^, the contribution from this 
extrapolation is very small and well defined. 

Table [T] lists the final results for the first moments measured at > 1 and evolved to a fixed 
(as indicated) for the more recent experiments. The quoted uncertainties include both statistical and 
systematic errors added in quadrature. The results from different experiments are in agreement within 
uncertainties, after taking Q^-evolution into account. The possible exception is the deuteron result of 
SMC [28] which is significantly lower than all others; this is due to a strongly negative contribution 
from the unmeasured low-x region which since has been ruled out by COMPASS data. In particular, 
the Bjorken sum rule is confirmed to within 7-15% of its value 0.182 ± 0.002 (evolved to the average 

= 5 GeV^ of the data listed in Table [I]); by far the largest remaining uncertainty stems from the 
low-x convergence of the integral. This agreement can be considered a successful test of Q^-evolution 
and the value of the strong coupling constant as- 

COMPASS and HERMES used their most recent result on the deuteron to extract an "experimental 
estimate" of the singlet axial charge ao (which equals AE in the MS scheme), yielding AE fa 0.35 ±0.06 
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for COMPASS (at = 3 GeV^) and AS ^ 0.33 ± 0.04 for HERMES (evaluated at = 5 QeV^). 
These numbers are somewhat higher than most recent NLO analyses (see Section 1?^ . but in agreement 
with similar analyses by other experiments. All results on Fi point towards a negative contribution 
to the integral from strange quarks and antiquarks of order ^ —0.1, while semi-inclusive results from 
HERMES I102j are consistent with As > in their measured range, see Sections I2.3|2.5[ 

3.3 The BC and ELT Sum Rules 

While the first moment of gi depends on and changes significantly from the DIS region to the real 
photon point, the Burkhardt-Cottingham sum rule for the first moment of g2 should apply at all values 
of > 0, as long as the elastic contribution is included in the integral. It can therefore be tested both 
by DIS data and at lower Q^, where a significant contribution to the integral comes from the nucleon 
resonance region. 

The SLAG E155 data discussed in Section [2^ yielded a first test of the sum rule for the proton and 
the deuteron. The data were integrated over the measured region, 0.02 < x < 0.8 at = 5 GeV^. The 
contribution from the unmeasured large-x region is negligible. For the low-x expansion, one can use the 
assumption that g2 follows the Wandzura-Wilczek form (Eq. [8OI) which yields a result independent of 
our knowledge of gi below x = 0.02. Under these assumptions, E155 found the integral for the proton 
to be —0.022 ± 0.008 and that for the deuteron as —0.002 ± 0.011, after averaging with the data from 
E143. While the proton result appears to be inconsistent with the BG sum rule at the 2.75a level, 
any firm conclusion depends strongly on the behavior of g2 at small x which is not known with enough 
precision and may not follow the Wandzura-Wilczek form. 

The same data were also used to estimate the value for the ELT integral (Eq. [201) using the approx- 
imation given in Eq. [HH with the neutron structure functions inferred from the measured proton and 
deuteron ones. The integral over the measured region is consistent with the expected value of zero, 
within errors of ±0.008. Here, the extrapolation to small x is less critical because of the extra factor of 
X in the integral. 



• JLabE94010 
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Figure 36: Results for the BG sum F2(Q^) for ^He |167j and the neutron [57]. The integrals 
over the measured region (filled circles) are compared with MAID model calculations. Light 
(dark) shaded bands indicate experimental (extrapolation) systematic errors. 

The most extensive measurement of the BG sum rule at smaller comes from an experiment using 
a longitudinally and transversely polarized '^He target in Hall A. The integral of Fj^^ = / g2^^{Q'^)dx 



50 



is plotted in the left panel of Fig. [36] for the measured region (from break-up threshold on up and cov- 
ering the resonance region; solid circles). For comparison, we also show a calculation using the MAID 
code |168j which is a unitary isobar model describing the nucleon resonance region, with parameters 
fit to meson photo- and electro-production data. Added to the MAID result is an estimate of the pd 
and ppn breakup contributions via a calculation of quasi-elastic single-proton knockout in the impulse 
approximation. The open circles show the experimental results after adding an estimated DIS contribu- 
tion. The solid squares correspond to the results obtained after adding the elastic contributions for ^He 
and are in good agreement with the expected value of zero within the systematic errors (indicated by 
the error bands for the experimental systematic error and the uncertainty on the low-x extrapolation). 

After applying some nuclear corrections, results for the neutron can be extracted from the ^He 
data. The result for T"!^ is plotted in the right panel of Fig. [36] in the measured region (solid circles), 
while the open circles include the elastic contribution and the open diamonds correspond to the results 
obtained after adding an estimated DIS contribution assuming g2 = for the neutron. Taking the 
difference between proton and deuteron data from El 55, one can estimate the result for the neutron 
at high (open square), which is consistent with zero but with a rather large error bar. The results 
from the Hall-A experiment for the neutron integral are the most precise data on the BC sum rule 
and are consistent with the expectation of zero, within systematic errors. The MAID parametrization 
agrees with the general trend but is slightly lower than the data in the measured (resonance) region, 
presumably because multi-pion final states are not included in the MAID code. 

The RSS experiment in Hall C [HOj also took data on g2 for the proton and the deuteron, at an 
average of about 1.3 GeV^. For the proton, they find a negative integral over the measured region 
(1 GeV < W < 2 GeV), but after adding the elastic contribution and an estimate for the small-a; 
part, they find a (preliminary) result for the BC sum very close to zero, with (small) statistical and 
(dominant) systematical errors of about ±0.005. This result is in disagreement with the SLAG result, 
but the unmeasured region (x < 0.29) is rather large and is significantly lower. Preliminary results 
on the deuteron are slightly negative but consistent with zero within an error of about ±0.015. 

3.4 Higher Moments and Higher Twist 

While the moments of the structure function gi provide information on the parton densities, the mo- 
ments of g2 are related to higher twist effects and provide information on quark-gluon interactions 
(correlations). The moment of the spin-structure functions given in Eq. [HHl 



is of special interest: at high it is a twist-3 matrix element and is related to the color polariz- 
abilities [169] , see Section 13.11 It can be calculated in lattice QGD |170j and a number of theoretical 
models [Wl [T50l [T5T] . 

Experimentally, due to the weighting, the contributions are dominated by the high-x region and 
the problem of low-x extrapolation is avoided. The SLAG E155 data [H] allowed the first extraction 
of ^2 and ^2 (excluding the elastic contribution). Gombining E155 results with earlier SLAG data, at 
an average of 5 GeV^ dl = 0.0032 ± 0.0017 and 4 = 0.0079 ± 0.0048. Gombining these data with 
the JLab g2 from Hall A E99-117 [62] a new value for the third moment ^2 was extracted at an average 

of 5 GeV^. Compared to the previously published SLAG result ^1], the uncertainty on ^2 has been 
improved by about a factor of 2 (see Fig. [571) . While a negative or near-zero value was predicted by 
Lattice QGD and most models, the new result for ^2 is positive. Also shown in Fig. [37] are the low 

(0.1-1 GeV^) results for ^2 from another Hall A experiment, E94-010 [57], which were compared 
to a Ghiral Perturbation Theory calculation |171] and a prediction based on MAID |172] . The elastic 
part, while negligible at high Q^, is significant at low {< 1 GeV^) and was subtracted from both 




(129) 
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Figure 37: Results for d^{Q^) from JLab Hall A [571 [62] and SLAG [H], together with Lattice 
QCD calculations |170j . 



theoretical calculations. The MAID calculation represents the low-Q^ data rather well, probably due 
to the x^-weighting which suppresses the contribution of higher mass final states. 

A new precision experiment to an average of 3 GeV^ is planned in Hall A in early 

2009 |173] . Further measurements |174j of (ig constant ranging from 3 to 5 GeV^ are planned at 
JLab after the 12 GeV energy upgrade. 

The Hall G RSS experiment [80] measured g2 on the proton and the deuteron and extracted (if = 
0.0072±0.0017 at a value of L3 GeV^. A more comprehensive measurement of g2 and ^2 is scheduled 
in Hall G |175] . It will cover a wide range from 2.5 to 6.5 GeV^. 

The higher-twist contributions to Fi can be obtained by a fit with an OPE series, Eq. ( 195|) . truncated 
to an order appropriate for the precision of the data. The goal is to determine the twist-4 matrix 
element f2- Once /i4 is obtained, /2 is extracted by subtracting the leading-twist contributions of 02 
and d2 following Eq. fllOOl) . To have an idea how the higher-twist terms (twist-6 and above) affect the 
twist-4 term extraction, it is necessary to study the convergence of the expansion and to choose the 
range in a region where fig term is not significant. This study is made possible only with the availability 
of the new low-Q^ data from JLab. 

Higher-twist analyses have been performed on the proton |176l 11771 1178j . the neutron [59] and the 
Bjorken sum {p — n) |179l llSOj . Fi at moderate was obtained from the JLab gi data with details 
described in Section 13. 5[ For consistency, the unmeasured low-x parts of the JLab and the world data 
on Fi were re-evaluated using the same prescription [59]. The elastic contribution, negligible above 

of 2 GeV^ but significant (especially for the proton) at lower values of Q^, was added using the 
parametrization of Ref . |181] . The leading-twist term fi2 was determined by fitting the data at > 5 
GeV^ assuming that higher twists in this region are negligible. Using the proton (neutron) data 
alone, and with input of a3(= gA) from the neutron beta decay and ag from hyperon decay (assuming 
SU(3) flavor symmetry), AE = 0.15 ± 0.07 for the proton analysis |176] and AE = 0.35 ± 0.08 for the 
neutron analysis [SH] were obtained. Note that there is a difference of nearly two standard deviations 
between these two results; this difference presumably comes from the different data sets used for the 
two analyses. 

The fit results using an expansion up to {1/Q^) in determining /i4 are summarized in Table [2j The 
extraction on p — n was performed twice: first |179j with the combined Hall A E94-010 neutron and 
Hall B EGIA proton data and then the second time [180] with the EGIB proton and deuteron data. In 
order to extract /2, shown in Table [3l the target-mass corrections 02 were evaluated using the Bliimlein- 
Bottcher NLO fit to the world data [133j for the proton and a JLab fit to the world neutron data, which 
includes the recent high precision neutron results at large x [611 [62]. The d2 values used are from SLAG 
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E155 gl] (proton) and JLab E99-117 [62] (neutron). 



Table 2: Results of /i4, /ig and /ig at = 1 GeV^ for proton, neutron and p — n. The 



uncertainties are first statistical then systematic. 



Target 


g2 (GeV2) 


/i4/M2 




/i8/M6 


proton 1176] 
neutron [59] 
p-n [1791 
p-n [180] 


0.6-10.0 
0.5-10.0 
0.5-10.0 
0.66-10.0 


-0.065±0.012± 0.048 
0.019±0.002± 0.024 
-0.060±0.045± 0.018 
-0.039±0.010± 0.026 


0.143±0.021 ±0.056 
-0.019±0.002± 0.017 
0.086±0.077± 0.032 
0.084±0.011± 0.024 


-0.026±0.008± 0.016 

0.00±0.00±0.03 
0.011±0.031 ±0.019 
0.047±0.026 



Table 3: Results of /2, Xe and xb at = 1 GeV^ for proton, neutron and p — n. The 



uncertainties are first statistical then systematic. 



Target 


/2 


Xe 


Xb 


p [176] 
n [59] 
p-n |179| 
p-n [180] 


-0.160 ±0.028 ±0.109 
0.034 ±0.005 ± 0.043 
-0.136 ±0.102 ±0.039 
-0.101 ±0.027 ±0.067 


-0.082 ±0.016 ±0.071 
0.031 ±0.005± 0.028 
-0.100 ±0.068 ±0.028 
-0.077 ±0.050 


0.056 ±0.008 ±0.036 
-0.003 ± 0.004 ± 0.014 
0.036 ±0.034 ±0.017 
0.024 ±0.028 



The fits were repeated varying the minimum threshold to study the convergence of the OPE 
series. The extracted quantities have large uncertainties (dominated by the systematic uncertainty) but 
are stable with respect to the minimal threshold when it is below 1 GeV^. The results do not vary 
significantly when the fig term is added, which justifies a posteriori the use of the truncated OPE series 
in the chosen range. In the proton case, the elastic contribution makes a significant contribution to 
the /ig term at low but this does not invalidate a priori the validity of the series since the elastic 
contribution affects mainly /ig and /xg remains small compared to yU4. We notice the alternation of 
signs between the coefficients. This leads to a partial suppression of the higher-twist effects and may 
be a reason for quark-hadron duality in the spin sector (see Section Hj). We also note that the sign 
alternation is opposite for the proton and neutron. Following Eq. (I103p . the electric and magnetic color 
polarizabilities were determined, see Table O We observe a sign change in the electric color polarizability 
between the proton and the neutron. We also expect a sign change in the color magnetic polarizability. 
However, with the large uncertainty and the small negative value of the neutron xb, if is difficult to 
confirm this expectation. 

Additional data from JLab experiment EGl are presently under analysis and will significantly im- 
prove the precision of the first moments and F^ around Q"^ = 1 GeV^, which is the crucial region for 
the extraction of higher-twist contributions. These data should therefore allow a more precise extraction 
of quantities like /2 and the color polarizabilities. 

3.5 GDH Sum Rule and Chiral Expansion 

In recent years the Gerasimov, Drell and Hearn (GDH) sum rule [1601 1161] at = has attracted 
much experimental and theoretical [BUI 1182] effort that has provided us with rich information. Experi- 
mental measurements of the GDH sum on the proton and deuteron were performed at Mainz [183] and 
Bonn [1841 1185] , and the contribution from the first resonance region and below was also measured at 
LEGS [186] . The results show that the GDH sum rule is satisfied for the proton, assuming a model [182] 
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for the unmeasured regions. A review on the overall status of the GDH sum rule for real photons has 
been published recently in this journal |187] . 

A generalized sum rule [188] connects the GDH sum rule with the Bjorken sum rule and provides a 
clean way to compare theories with experimental data over the entire range. The (generalized) GDH 
sum rule and similar spin sum rules relate the moments of the spin-structure functions to the nucleon's 
static properties (e.g., its anomalous magnetic moment) and to real or virtual Compton amplitudes 
(see Eqs. I116tll22tfl25tll27l) . which can be calculated theoretically. Several papers [60l I164[ 1182] provide 
comprehensive reviews on this subject. 




Figure 38: Results of the GDH sum I{Q^) for ^He |167j (left) and for the neutron [56] 
(right). The ^He GDH results are compared with the MAID model plus a contribution 
from quasi-elastic nucleon knock-out. The neutron GDH results are compared with xPT 
calculations of Ji |171] (dashed line) and Bernard et al. [189] (shaded area). The MAID 
model calculation [172] is represented by a solid line. Data from HERMES [19Uj are also 
shown. 

Fig. [38] shows the extended GDH integrals I{Q^) = JI^[ai/2{Q^) - <y-i/2{Q'^)]dv /v for ^He (left) and 
for the neutron (right), which were extracted from Hall A experiment E94-010 [1671 from break-up 
threshold for ^He and from pion threshold for the neutron io W = 2 GeV. The uncertainties, when 
visible, represent statistics only; the systematics are shown by the shaded bands. The solid squares 
include an estimate of the unmeasured high-energy part. The corresponding uncertainty is included in 
the systematic uncertainty band. 

The ^He results rise with decreasing Q^. Since the GDH sum rule at = predicts a large 
negative value, a drastic turn around should happen at lower than 0.1 GeV^. A simple model 
using MAID [172] plus quasi-elastic contributions estimated from a PWIA model [191] indeed shows 
the expected turn around. The data at low should be a good testing ground for few-body Chiral 
Perturbation Theory calculations. 

The neutron results indicate a smooth variation of I{Q^) to increasingly negative values as varies 
from 0.9 GeV^ towards zero. The data are more negative than the MAID model calculation [T72], again 
presumably due to the contribution from multi-pion final states. (Since the calculation only includes 
contributions to /(Q^) for W < 2 GeV, it should be compared with the open circles.) The GDH sum 
rule prediction, J(0) = —232.8 /ib, is indicated in Fig. [381 along with extensions to > using two 
next-to-leading order xPT calculations, one using the Heavy Baryon approximation (HB^PT) [171] 
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Figure 39: Results on Ti{Q'^) for p, d, n, and p-n from JLab Hall A [57] and CLAS EGla [71] 
and EGlb [I92], together with data from SLAG and HERMES. The slopes at = g 
predicted by the GDH sum rule are given by the red dotted lines. The solid (short-dashed) 
lines are the predictions from the Burkert-Ioffe (Soffer-Teryaev) parametrizations |194[ 1193] . 
The leading twist Q^-evolution of the p-n moment is given by the grey band. The insets 
show comparisons with xPT calculations. The dot-dashed lines (green lines and bands) at 
low Q"^ are the next-to- leading order xPT predictions by Ji et al. |171j (Bernard et al. |189j ). 



(dashed line) and the other based on Relativistic Baryon xPT (RB^PT) |189] . The shaded band shows 
the RB^FT calculation including resonance effects [189j . which have an associated large uncertainty 
due to the resonance parameters used. 

In Fig.[39]we present the existing data set on the related quantity Ti^Q"^) (Eq. ll20| ) at low to moderate 
Q^. The upper two panels show the data on the proton and the deuteron from EGla [71] together with 
new results from the Hall B EGlb experiment [192] . The bottom row shows published results for the 
neutron from Hall A [57] and for the proton-neutron difference obtained by combining all of these data 
sets. Also shown are data from SLAG E143 [3l] and E155 [36] as well as from HERMES [S]. The 
inner error bars indicate the statistical uncertainties while the outer ones are the quadratic sums of the 
statistical and systematic uncertainties, except for the lower-right panel where the systematic errors are 
indicated by a shaded band. 
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As — 0, the slopes of Fi (shown as dotted hnes in Fig. [39]) and of Ttt are the same, both given by 
the GDH sum rule (see Eqs. I106mi0lll20|12ip . The behavior at low can be calculated with xPT. We 
show results of calculations by Ji et al. |171] using HBxPT and by Bernard et al. |189] without and with 
the inclusion of vector mesons and A degrees of freedom. The calculations are in reasonable agreement 
with the data at the lowest values of 0.05 - 0.1 GeV^. At moderate and large data are compared 
with two phenomenological model calculations |193[ 1194] . The model by Soffer and Teryaev |193] 
assumes a smooth variation with of the integral / [gi{x) + g2{xy\dx (which is constrained by DIS 
data at large and by sum rules at the photon point) and subtract the contribution from f 2 (using 
the EC sum rule) to extract TiiQ"^). The curve by Burkert and loffe |194] uses a parametrization of the 
resonance region (similar to MAID) and adds a Vector Meson Dominance-inspired ansatz to connect 
the GDH sum rule at the photon point with the DIS data for Fi. Both models describe the data well. 

The lower-right panel in Figure [39] shows the moment of — (7", the Bjorken integral. This is of 
special interest because it contains contributions only from the flavor non-singlet (or isovector) part, 
and the contribution from the A(1232) resonance cancels |195] . The data at high were used to test 
the Bjorken sum rule as described in Section [3l2l They were also used to extract a value of the strong 
coupling constant, a^, assuming the validity of the sum rule. The leading-twist pQCD evolution of the 
Bjorken sum rule is shown by the grey band. It tracks the data down to surprisingly low Q^, which 
indicates an overall suppression of higher-twist effects (see Section [374]) . The new JLab data at low 
provide interesting information in the low energy region, where the strong interaction is truly strong 
and non-perturbative. 
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Figure 40: Results for a^^g^ |196] extracted from the generalized Bjorken sum, together with 
other world data on effective strong couplings. 



A new attempt [196] 1197] was made to extract an effective strong coupling, in the low region 
(Figure [10]) • The extracted which by definition must converge towards a finite value as — > 0, 
shows a clear trend of weakening Q^-dependence with decreasing Q^. With the GDH sum rule as a 
constraint at = 0, a. model fit to the extracted ag^g^ shows a small Q^-slope at the origin. This is 
consistent with a conformal behavior, which may be important for any attempt to apply AdS/CFT |198] 
for the strong interaction in the low-energy region. 

The generalized spin polarizabilities provide additional benchmark tests of xPT calculations, and 
their measurement is an important step in understanding the dynamics of QCD in the low region. 
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Since they have an extra \jv^ weighting compared to the first moments, these integrals have less 
contributions from the large-i^ region and converge much faster, which minimizes the uncertainty due 
to the unmeasured region at large v. 

Generalized spin polarizabilities have been evaluated with next-to-leading order xPT calculations [199[ 
I189j . One issue in the xPT calculations is how to properly include the nucleon resonance contributions, 
especially the A(1232) resonance. As was pointed out in [199[ 1189] , while 70 is sensitive to resonances, 
8i,T is insensitive to the A resonance. 




Figure 41: Results for the neutron spin polarizabilities 70 (top panel) and bi,T (bottom 
panel). Solid squares represent the results with statistical uncertainties. The light bands 
represent the systematic uncertainties. The dashed curves represent the heavy baryon xPT 
calculation |199] . The dot-dashed curves and the dark bands represent the relativistic baryon 
xPT calculation without and with |189] the A and vector meson contributions, respectively. 
Solid curves represent the MAID model [172] . 



The first results for the neutron generalized forward spin polarizabilities 7o(Q^) and 5lt(Q^) were 
obtained at Jefferson Lab Hall A [5H]. The results for 7o (Q^) are shown in the top panel of Fig. SH 
The statistical uncertainties are smaller than the size of the symbols. The data are compared with a 
next-to- leading order (0(p^)) HBxPT calculation [199] , a next-to- leading order RB^PT calculation and 
the same calculation explicitly including both the A resonance and vector meson contributions [189j . 
Predictions from the MAID model [172] are also shown. At the lowest point, the RB^PT calculation 
including the resonance contributions is in good agreement with the experimental result. For the HB^PT 
calculation without explicit resonance contributions, discrepancies are large even at = 0.1 GeV^. 
This might indicate the significance of the resonance contributions or a problem with the heavy baryon 
approximation at this Q^. The higher data point is in good agreement with the MAID prediction, 
but the lowest data point at = 0.1 GeV^ is significantly lower. 

Results for 70 on the proton has been recently submitted for publication [192] . They show significant 
disagreement with both xPT calculations [199[ 1189] . An isospin separation of 70 was performed and 
discussed in Ref. |180j . The isoscalar (70^") and isovector (7o~") combinations also disagree with xPT, 
in spite of the fact that for the latter the contribution from the A(1232) resonance cancels. 
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Since 6lt is insensitive to this resonance contribution for each nucleon individually, it was believed 
that 6lt should be more suitable than 70 to serve as a testing ground for the chiral dynamics of 
QCD [199^ I189j . The bottom panel of Fig. HT] shows 6lt compared to xPT calculations and the MAID 
predictions. While the MAID predictions are in good agreement with the results, it is surprising to see 
that the data are in significant disagreement with the xPT calculations even at the lowest (5^ 0.1 GeV^. 
This disagreement ("5lt puzzle") presents a significant challenge to the present Chiral Perturbation 
Theory. 

New experimental data have been taken at very low Q^, down to 0.02 GeV^ for the neutron (^He) [53] 
for both longitudinal and transverse target polarizations, but only for longitudinal target polarization for 
the proton and the deuteron |200] . Preliminary results just became available for the neutron. Analysis 
is underway for the proton and deuteron data. These results will provide benchmark tests of xPT 
calculations at the kinematics where they are expected to work. A new proposal [201j was recently 
approved to measure g2 with a transversely polarized proton target in the low region. It will provide 
an isospin separation of the spin polarizabihties to shed light on the "5^^" puzzle. 

4 Quark-Hadron Duality in Spin Structure 

In the previous sections, we have interpreted spin structure function data and their integrals either in 
terms of quark and gluon degrees of freedom (Sections l2.1H2.6[l3.1H3.4() or in terms of low-energy effective 
theories (Sections 12. 7[ 13.51) . While quarks and gluons are the fundamental degrees of freedom for the 
theory of strong interactions, QCD, a description in terms of their hadronic composites (nucleon, nucleon 
resonances and light mesons) is often more economical at lower resolution Q^, where perturbative QCD 
is not applicable and calculations are extremely difficult. At intermediate resolution, models like the 
constituent quark model may be able to describe some of the observed phenomena. 

Both from a theoretical and a practical point of view, it is interesting to investigate where and how 
these different pictures of the nucleon connect and overlap. In particular, if we can find a kinematic 
region where both the quark-parton description and the hadronic description are found to be valid to 
some degree, we can potentially expand the kinematic range of experiments that attempt to extract 
information on quarks in the nucleon. At the same time, we could gain new insights on the transition 
from quasi-free quarks at short length scales to their confinement at larger distances. 

Such "dual" descriptions of experimental data have been successful in cases including e+e~ annihi- 
lation, semi-leptonic heavy meson decays and lepton scattering; for a review see |202j . In particular, 
in the 1970's Bloom and Gilman |203] found that the unpolarized structure function F2{x,Q'^) in the 
nucleon resonance region {W < 2 GeV) at moderate Q^, when averaged over suitable intervals in x, 
was quite close to the deep inelastic, scaling structure function ^2(0;) measured at much higher 
and W but the same average x. This finding has been confirmed and expanded in much greater de- 
tail by measurements at Jefferson Lab [204j . The agreement between resonance region and high-Q^ 
data is typically improved if one accounts for target mass effects and the logarithmic Q^-evolution of 
deep inelastic structure functions. Nachtmann [205] has shown how to handle the finite mass of the 
nucleon target exactly; an approximate method simply employs the modified Nachtmann scaling vari- 
able ^ = 2a;/ (1 + a/1 + 4M^x^/(5^)- This has been extended to the polarized case by Piccione and 
Ridolfi [T28] . 

In the framework of the OPE, a quark-parton description of structure function moments at interme- 
diate should in principle be always possible if one includes contributions of sufficiently high twists. 
This so-called global duality then becomes simply the observation that these higher-twist contributions 
are either small or cancel in integrated observables [206j . Experimentally, this seems to be the case for 
unpolarized structure functions of the proton; however, it is non-trivial whether the same is true for 
polarized structure functions and for neutron structure functions. In particular, the polarized struc- 
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ture function gi might not necessarily exhibit local duahty (for a restricted range in x), since at low 
and W it is dominated by the A(1232) resonance, which should have a negative spin asymmetry 
Ai ^ —0.5 (see Section [T7|) . while at large and x one expects this asymmetry to approach unity 
(see Section [2^ . Therefore, a detailed study of the limits of validity of duality for polarized structure 
functions is of great interest. 

First measurements of spin structure functions in the resonance region [20l [39] found that, with the 
exception of the A(1232) resonance, duality seemed to hold approximately within errors, especially for 
the higher data of E143 [39]. Early data from the Jefferson Lab EGl program with CLAS [701 [7T] 
confirmed that local duality does not work well at low Q^, with significantly higher statistics. The 
paper by Yun et al. (TDj explicitly looked at the approach of the structure function gf{^, Q"^) towards 
the scaling limit for increasing and found that this limit had not been reached yet up to = 1 
GeV^. The HERMES collaboration [207j averaged their measured asymmetry A\ over the resonance 
region 1 GeV < W < 2 GeV for 5 different bins in x and and found the results in reasonable 
agreement with DIS results, within their relatively large errors. 
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Figure 42: Data on the spin structure functions gi{x,Q^) of the proton and the deuteron 
from Jefferson Lab's Hall B (left panel) and its average over the resonance region (right 
panel). Prominent resonances are indicated by arrows in the left panel. The hatched curves 
represent the range of extrapolated DIS results from modern NLO fits (GRSV and AAC, 
see Section 123]) , evolved to the of the data and corrected for target mass effects. The 
open circles in the right panel include the elastic contribution, while the filled circles are 
only integrated over W > 1.08 GeV. 

By far the most detailed study of duality in the spin structure function gi has been published by the 
EGl collaboration [73]. As shown in Fig. HJ] one observes a clear trend of strong, resonant deviations 
from the scaling curve at lower Q^, towards a pretty good agreement at intermediate Q^. The integral of 
gi over the whole resonance region begins to agree with the NLO results above ~ 1.7 GeV^. At lower 
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Q^, that integral rises above the NLO extrapolation if one includes the elastic contribution. This shows 
that the elastic peak somewhat overcompensates the negative contribution from the A(1232) resonance; 
however, leaving out the elastic part of the integral leads to integral values much below the NLO 
curves. The results on the proton and deuteron from EGl [73] thus indicate a much slower approach to 
"global" duality for the polarized structure function gi than has been observed for unpolarized structure 
functions. 

A more detailed examination shows that gi averaged over the first resonance region (1.08 < W < 1.38 
GeV), which includes the A(1232) resonance, is substantially below the NLO expectations even to the 
highest measured by EGl (5 GeV^). On the other hand, the three resonances (Pii(1440), Sii(1535), 
Di3(1520)) in the second resonance region (1.38 < W < 1.58 GeV) all have dominant A1/2 transition 
amplitudes already at rather modest Q^, so that the averaged gi in this region tends to overshoot the 
NLO limit. The cancellation of these two opposite effects, together with a rather quick approach to 
"local duality" for the higher resonance region (W > 1.56 GeV) leads to the approach to global duality 
seen in Fig. | 
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Figure 43: New data on the spin structure functions gi{x,Q'^) of '^He in the resonance 
region from Jefferson Lab's Hall A (solid circles, left panel), together with NLO curves for a 
combination of proton and neutron gi with proper polarization factors to mimic the nuclear 
spin structure function and evolved to the of the data. The right panel show these 
data (solid squares) averaged over the resonance region for ^He (top half) and corrected for 
nuclear effects to get the partial integrals for the neutron (bottom half). The NLO curves 
shown represent the range of extrapolated DIS results from modern NLO fits (Bliimlein and 
Bottcher, GRSV, LSS and AAC, see Section [2751) . evolved to the of the data and corrected 
for target mass effects. 



The data taken by the RSS collaboration in Hall C (see Section [277|) corroborate these observations 
and add more precise data points for ~ 1.3 GeV^. They find that the integral over their measured 
region, 1.09 GeV < W < 1.91 GeV, is about 17% ± 8% lower than the same integral over NLO 
evaluations of g^^ evolved to = 1.3 GeV^, even after applying target mass corrections. 

Very recently, new results have become available from a duality study in Hall A, using a "^He tar- 
get [M]- By combining NLO fits for the spin structure functions g\ and g^ with the proper polarization 
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factors, one can directly compare these fits with the data on the nucleus ^He (left panel of Fig. H3|) . Once 
again, at the lowest point the data oscillate around the NLO curves and show a strong deviation in 
the region of the A(1232) resonance at least out to = 1.8 GeV^. At higher Q^, both the data and 
the NLO fits are very close to zero, in contrast to g\ which is markedly positive for the extrapolated 
NLO fits. The partial integrals (up to = 1.91 GeV) of these data are shown in the right panel of 
Fig. B3|) . Except for the lowest point, good agreement with the integrals of NLO fits (both for ^He 
and directly for the neutron) is seen, confirming the onset of duality around = 1.7 GeV^ for both 
protons and neutrons. 

5 Summary and Outlook 

After 30 years of dedicated experiments, we are beginning to accumulate a rather detailed picture of 
a fundamental aspect of the nucleon, namely its spin structure. In spite of early worries (after the 
initial "spin crisis"), it appears that quark helicities do contribute a small but significant fraction of the 
(longitudinal) spin of the nucleon, up to one third. Most importantly, the fundamental Bjorken sum 
rule [S] and its perturbative QCD evolution seem to be at least consistent with the data. In fact, none 
of the existing experiments show any features that would contradict pQCD in its realm of applicability. 

Most pQCD-based analyses of the data agree fairly well on the contribution of various quark flavors 
to the proton spin, although there is still some controversy on the role played by strange sea quarks. 
We also don't know yet whether polarized up and down anti-quark densities show the same difference 
as the unpolarized ones. On the other hand, new high-precision data from Jefferson Lab constrain 
very nicely the polarized densities of valence up and down quarks at higher x. These data indicate 
that relativistic quark models, with one-gluon exchange corrections |lU9j , and augmented by pion cloud 
contributions |2U8j . yield a rather good description in the valence region where sea-quarks contribute 
little. The emergent picture (also from recent NLO analyses, see Section 12.51) is one where valence 
quarks (Agy = Aq — Aq) carry roughly the expected fraction (2^ 60%) of the nucleon spin, while the (on 
average) negative helicity of sea quarks reduces this to about 30-35%. Future experiments |209U210[I211] 
at Jefferson Lab, after the approved energy upgrade to 11-12 GeV, will extend our knowledge of polarized 
quark densities out to x > 0.8 and will decisively test predictions from pQCD and QCD-inspired models. 
These same experiments, as well as presently scheduled 6 GeV ones |212[ 1213] , will further improve the 
precision with which we can study scaling violations and therefore extract polarized gluon and sea- 
quark densities. Additional information on individual quark and antiquark flavor contributions to the 
nucleon spin will come from future experiments at RHIC (in particular from direct production 
at 500 GeV center-of-mass energy) and the new FAIR facility in Darmstadt, Germany |214] as well as 
semi-inclusive measurements with COMPASS and at Jefferson Lab. Much additional information on the 
small-x behavior of quark- and antiquark polarization densities could be gotten from the electron-ion 
collider (EIC) that has recently been proposed |215j . see below. 

In addition to logarithmic violations of scaling, expected from pQCD, we also see some evidence for 
non-zero twist-3 and twist-4 matrix elements; however, the overall higher twist corrections to moments 
of spin structure functions seem rather modest even at fairly low Q^, indicating a partial cancellation 
of contributions from different orders in the twist expansion. At these lower values of Q^, we see a 
transition from partonic (quark and gluon) degrees of freedom to hadronic ones, without any dramatic 
break in the rather smooth transition region. Thanks to the vast data set in the nucleon resonance 
region, we can better constrain transition amplitudes and test this quark-hadron duality in detail. 
Averages over the resonance region seem to agree fairly well with extrapolations of deep inelastic data 
over the same range in x, albeit not to quite as low as in the unpolarized case. Future measurements 
at Jefferson Lab with 6 GeV [1751 1173[ 1216] and 12 GeV |174j beams will add to our knowledge of 
higher-twist matrix elements, and of the second spin structure function ^2 and its moments. 
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At the lower end of the scale, the Gerasimov-Drell-Hearn sum rule at the real photon point 
seems to be well confirmed. However, it is too early to state with certainty whether the effective theory 
of QCD at low energies, xPT, yields a good description of the data very close to that point. Some 
tentative successes in describing the first moment Fi of gi will be tested very soon, when the new data 
on the neutron (^He) , the proton and the deuteron collected in Halls A and B at Jefferson Lab become 
available [651 1200j . For the polarizabilities like 70 and 5lti which are given by higher moments of (71^2, 
the agreement at present is unsatisfactory; it may be necessary to go to a higher order in xPT before 
agreement is reached. New data on these quantities will also be collected at Jefferson Lab |201j . 




Figure 44: gi{x,Q'^) at various values of attainable at EIC for positive and negative AG |217] . 

One of the most important remaining open questions is the provenance of the remaining 2/3 of the 
nucleon spin that is not carried by quark helicities. We know now that the gluon contribution cannot 
be very large, but its precise magnitude and shape are not known yet, and it could still be an important 
fraction of the total. Improved statistics from the direct measurements at RHIC and COMPASS will 
help, especially with expanded kinematics (higher and lower beam energies, different rapidity and pt 
ranges, and new final state channels) which can begin to separate the low-x and high-x behavior of AG. 

Additional information will come from future experiments of the spin structure functions at both 
extremes of the scale, including the upcoming Jefferson Lab experiments with the future 11 GeV 
beam |209j . and further COMPASS data. In this context, it is an interesting question whether future 
DIS measurements alone can distinguish between positive and negative AG(a;). They can, indeed, but 
this requires a very high energy lepton- nucleon collider like the EIC proposal |215] . This can be seen 
in Fig. HH which shows g\{x) at very small x for 1 < < l,OOOGe^^ for the two signs of AG(x). 
Clearly the behavior of gi{x) at small x and large is quite different in the two cases. 

A collider of the EIC type would also have a dramatic effect in reducing the uncertainties in the 
polarized parton densities. This is illustrated in Figs. |45|^ where it can be seen that the improvement 
even at moderate to large x, especially for As and AG, is remarkable. 

Finally, the contribution from quark angular momentum is also an important ingredient in the 
total spin balance of the nucleon. While a direct measurement is not available, one can learn much 
about the transverse distribution and motion of quarks from semi-inclusive measurements of single 
spin asymmetries. This is a fairly new field, with a very rich potential and a rapidly growing body of 
experimental data, but lies outside the scope of this article. It will be addressed by an upcoming review 
in this journal. 

A second avenue towards a full accounting of the nucleon spin lies in measurements that are sen- 
sitive to Generalized Parton Distributions (GPDs), in particular Deeply Virtual Compton Scattering 
(DVCS). Moments of certain combinations of GPDs can be related to the total angular momentum 
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Figure 46: Impact of two versions of EIC on the s and G uncertainties |217] . 



(spin and orbital) carried by various quark flavors, as expressed in Ji's sum rule [H I218j . Again, semi- 
nal experiments in this area have already collected data that indicate that the necessary assumption of 
factorization (and the validity of the handbag diagram) seems to be fulfilled, and even give some model- 
dependent indication of the value of u- and d-quark angular momentum |219] . Future measurements at 
COMPASS and an extensive program at the energy-upgraded Jefferson Lab will add substantially to 
this database. Ultimately, an electron-ion collider could map out GPDs in detail at significantly lower 
X than fixed-target experiments. 

In summary, the interest in the nucleon spin structure continues unabated, with many novel ex- 
perimental and theoretical approaches towards the ultimate understanding of all contributions to the 
nucleon's spin. An extensive and rich experimental program lies already ahead of us at existing facilities 
like CERN (COMPASS), RHIC and Jefferson Lab (both during the remaining years of 6 GeV running 
and in the new era of 11-12 GeV beams). Ultimately, this field would benefit tremendously from a new 
collider of both polarized ions and polarized leptons, which is now in the planning stage. 
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